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0\ ■ ABSTRACT 

^: 

, ^ \ We have obtained moderate-resolution {R = 800-1200) i^T-band spectra for ~ 100 

stars within and surrounding the cloud core of p Oph. We have measured spectral 
types and continuum veilings and have combined this information with results from 
new deep imaging. Using the latest evolutionary tracks of D'Antona &; Mazzitelli to 
interpret the H-R diagram for p Oph, we infer ages ranging between 0.1 and 1 Myr for 
the Class II and III sources (i.e., those that have emerged from their natal cocoons). A 
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OO . few stars may be slightly older. The IMF peaks at about 0.4 Mq and slowly declines 

\ to the hydrogen burning limit with a slope of ~ —0.5 in logarithmic units (Salpeter is 

\ -1-1.35). Our lower limits on the numbers of substellar objects demonstrate that the 

' IMF probably does not fall more steeply below the hydrogen burning limit, at least 

' down to ~ 0.02 Mq. The derived IMF is consistent with previous findings that the 

Q^' p Oph IMF is roughly flat from 0.05 to 1 Mq. The exact shape of the mass function 

Q . remains a function of the theoretical evolutionary tracks and, at the lowest masses, the 

\ We then make the first comparison of mass functions of stars and pre-stellar clumps 

(Motte, Andre, & Neri) measured in the same region. The similar behavior of the two 
^ mass functions in p Oph supports the suggestion of Motte et al. and Testi & Sargent 

;h ■ that the stellar mass function in young clusters is a direct product of the process of 

^ cloud fragmentation. 

We have also studied the very young and often still embedded Class I and flat 
spectrum objects. After considering the effect of extinction on the SED classifications 
of the sample, we find that ~ 17% of the p Oph stars are Class I, implying ~ 0.1 Myr for 
the lifetime of this stage. In spectra separated by two years, we observe simultaneous 
variability in the Br7 emission and ET-band continuum veiling for two stars, where the 
hydrogen emission is brighter in the more heavily veiled data. This behavior indicates 
that the disk may contribute significantly to continuous K-hand emission, in contrast 
to the proposal that the infalling envelope always dominates. Our detection of strong 
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2 /iiri veiling [vk = 1-4) in several Class II and III stars, which should have disks but 
little envelope material, further supports this proposition. We also detect absorption 
features in the spectra of ~ 25% of Class I and flat-spectrum sources, demonstrating 
the feasibility of studying the photospheres of extremely young protostars. 

Subject headings: infrared: stars — stars: evolution — stars: formation — stars: 
low-mass, brown dwarfs — stars: luminosity function, mass function — stars: pre-main 
sequence 

1. Introduction 

The birth of low-mass stars and brown dwarfs as well as the evolution of circumstellar disks 
and planetary systems can best be studied in the nearest {d < 500 pc) young (age < 10 Myr) 
stellar clusters and associations. For example, the distribution of masses in these populations 
directly reflects the initial mass function (IMF). In contrast, in regions like the local field or 
globular clusters, the present day mass function is the product of a complex evolution. 

Until recently, these advantages of young clusters were difficult to exploit because of 
extinction. However, by complementing optical techniques with modern infrared (IR) imaging and 
spectroscopy, we can estimate luminosities and spectral types for complete, well-defined samples 
of young objects within obscured star forming regions. Due to the luminous nature of newborn 
substellar objects, such data can be converted to mass functions that reach below the hydrogen 
burning limit. 

The p Ophiuchi dark cloud core contains a very young (< 1 Myr), nearby {d = 160 pc), 
compact {D ~ 20') population of ~ 100 low-mass stars, making it ideal for IMF studies. Due 
to the large extinction within the cloud, optical spectral types and photometry (Bouvier & 
Appenzeller 1992; Martm et al. 1998) and soft X-ray data (Montmerlc et al. 1983; Casanova et 
al. 1995) are available for only a small fraction of the stars. Until recently, the p Oph members 
were studied primarily through their near- to mid-IR spectral energy distributions (SEDs), which 
were observed with single element detectors on the ground and aboard the IRAS satellite (Elias 
1978; Lada & Wilking 1984; Wilking, Lada, & Young 1989; Greene et al. 1994). Following the 
initial observations of p Oph with near-IR arrays (Barsony et al. 1989; Rieke, Ashok, Sz Boyle 
1989), Greene & Young (1992) surveyed 650 arcmin^ of the cloud at J, H, and K, providing a 
relatively sensitive (K ~ 13) photometric census of the stellar population. To reach young objects 
below the hydrogen burning limit, Rieke & Rieke (1990) obtained deep images {K ~ 15.5) of 
200 arcmin^ within the cloud core {Ay ^ 50, Wilking & Lada 1983) and identified several brown 
dwarf candidates. The substellar nature of these sources has been supported by ISO mid-IR 
photometry (Comeron et al. 1998, hereafter CRCTL) and confirmed through spectroscopy in the 
optical (Luhman, Liebert, & Rieke 1997) and IR (Williams et al. 1995; Wilking, Greene, &; Meyer 
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1998, hereafter WGM). Comeron et al. (1993) (hereafter CRBR) supplemented the data of Rieke 
Sz Rieke (1990) with JHKL' aperture photometry of a complete sample of sources projected 
against the dense cloud core. They introduced a method of fitting the data to estimate individual 
source properties and deduced a logarithmically flat IMF (the Salpeter index is +1.35) down to 
and below the hydrogen burning limit. Strom, Kepner, & Strom (1995) used an alternate approach 
introduced by Meyer (1995) that also indicated a flat IMF. Williams et al. (1995) showed that 
the techniques of CRBR and Strom et al. (1995) agree at the detailed level of individual source 
properties as well as overall IMF slope. 

Following these photometric studies, near-IR spectroscopy has been applied in studying other 
properties of the embedded stellar population in p Oph. After measuring IC-band spectral types 
for 19 stars {K < 10.5, R = 500), Greene &: Meyer (1995) constructed a Hertzsprung-Russell (H-R) 
diagram in which the theoretical evolutionary tracks of D'Antona &l Mazzitelli (1994) (hereafter 
DM94) implied an age of ~ 0.5 Myr for the cluster. Greene k. Lada (1996) also observed ~ 100 
young stars in Taurus and p Oph [R = 500) with broader wavelength coverage (1.15-2.42 /xm). 
The appearance of the spectra was closely correlated with SED class, indicating, as expected, 
that reddening increases and photospheric features weaken from Class III to Class I sources 
(Lada 1987) and thus confirming the previous conclusions drawn from spectroscopy of 2.3 pm 
CO absorption in young p Oph members (Casali &; Matthews 1992). This trend is consistent 
with the large obscuration and strong continuum veiling expected from circumstellar material in 
the presumably younger, earlier SED classes. Finally, Casali &; Eiroa (1996) and Greene & Lada 
(1997) used high resolution spectra {R = 17,000 and 21,000, respectively) to measure continuum 
veilings and line widths for a small sample of flat-spectrum sources in p Oph and other young 
clusters. They demonstrated that the line profiles arc consistent with absorption originating in 
very young stellar photospheres rather than circumstellar disks, an important result that supports 
the validity of our classification of young stellar populations through IR spectra. 

In this paper, we expand on the previous IR spectral classification in p Oph and continue a 
program started in L1495E (Luhman &; Rieke 1998, hereafter LR) and IC 348 (Luhman et al. 
1998b, hereafter LRLL). We have obtained K-band spectroscopy of ~ 100 sources, combining 
a magnitude-limited sample in the cloud core {K ^ 12) with a representative population from 
the outer regions of the cluster (K <^ 11). The moderate spectral resolution (R = 800-1200) 
of the data allows both reasonably accurate spectral types (±1-2 subclasses) and sensitivity to 
low-mass members of p Oph (> 0.1 Mq). With the measured spectral types we construct an 
H-R diagram and use theoretical evolutionary tracks to derive an IMF and star formation history 
for the emerged members of the cluster (predominantly Class II and III). Likely substellar cloud 
members examined by Luhman et al. (1997), WGM, and CRCTL are added to the mass function, 
extending the IMF to lower masses. We then use a new K-hand luminosity function (KLF) of the 
cloud core to develop a model of the reddened background population and identify likely cluster 
members lacking spectroscopy. After accounting for these sources, we arrive at an IMF complete 
to ~ 0.08 Mq, which is compared to results from other young clusters. 
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Using the mid-IR photometry of previous authors, we revise the SED classifications of the 
p Oph population by dereddening the near- to mid -IR spectral indices with guidance from our 
near-IR spectra. With the new classifications, the relative numbers of sources in the different 
categories are consistent with expectations that the dust embedded stage should have a very short 
duration. The broad-band photometric data combined with the measurements of the K-haiid 
continuum veiling and Br7 emission present new constraints on the nature of the IR excess 
emission in young stars. In particular, simultaneous variability of Br7 emission and veiling is seen 
in two stars, possibly indicating the importance of disks as sources of 2 //m excess emission. In 
addition, we find that photospheric absorption features are detectable in a number of Class I and 
flat-spectrum sources, while significant veiling is measured in several Class II objects. Both the 
placements on the H-R diagram and the distributions of SED classes indicate a slight increase in 
average age between the cloud core and surrounding region. 

2. Observations 

We performed i^-band spectroscopy on sources in p Oph using the near-IR long-slit 
spectrometer FSpec (Williams et al. 1993) at the Steward 2.3 m Bok Reflector on Kitt Peak on 
1994 July 15 and 1995 May 16, 17, 19 and at the Multiple Mirror Telescope (MMT) on Mount 
Hopkins on 1994 July 2-3, 1996 May 28-31, and 1996 June 3-6. Most of the spectra were obtained 
in 1996 with a grating providing a two-pixel resolution oi R = A/ A A = 1200, while the data from 
1994 and 1995 have a resolution of i? = 800. The observations and data reduction procedures were 
identical to those described by LRLL. We selected for spectroscopy ~ 30 sources with K ^ 12 
appearing in IR images of the cloud core by CRBR. With the addition of spectra from Williams 
et al. (1995), our sample contains ~ 40 sources and is nearly complete to K = 12. In the region 
surrounding the core and within Figure |l|, we observed ~ 70 sources with K ^ 11 found in Greene 
& Young (1992). A few bright stars outside of the region in Figure |l] were also observed, taken 
from Wilking, Lada, & Young (1989) and Greene et al. (1994). 

The spatial distribution of the spectroscopic sample is illustrated in Figure |l| and the spectra 
are presented in Figs. Line identifications were made by comparison to the high-resolution 
(i? > 45, 000) spectra of cool stars obtained by Wallace & Hinkle (1996). 

Imaging at Ks (1.9-2.3 pm) was performed at the Bok Refiector on the nights of 1998 April 
12-15 using the Steward Observatory NICM0S3 256 x 256 near-IR camera at a plate scale of 
0'.'64 pixel"^, corresponding to a total field of 2'.7 x 2'.7. In 43 pointings, we observed the same 
region covered by CRBR (outlined by the dashed lines in Figure |l]), in addition to a few extra 
positions, for a total of ~ 250 arcmin^. Followup observations were made with standard J and H 
during 1998 May 8-10 and July 14. We selected pointings that contained several faint, possibly 
low-luminosity i^T-band sources (outlined by dotted lines in Figure Q). At each position, we 
obtained images with dithers of 5" in a 4 x 4 grid, facilitating efficient fiat-fielding and removal of 
bad pixels. Each exposure consisted of two coadded 30 s frames, producing a total exposure time 
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of 16 minutes. A flat field was constructed by median combining frames that lacked bright sources. 
For cahbration, we observed a photometric standard star, SR3, periodically during the night. 
These exposures were taken out of focus to guard against saturation. After dark-subtracting and 
flat-fielding the images, the 16 frames at a given position were shifted and combined into one 
image. The tasks DAOFIND and PHOT under the package APPHOT were used in measuring 
the stellar coordinates and extracting photometry. The density of stars was low enough that we 
could measure the background emission in an annulus around each star and subtract it from the 
photometry. The data were calibrated with the CIT photometry of SR3 from Elias et al. (1982) 
and are essentially in the CIT system. The K images are complete to -fC ~ 17, based both on 
the detection of many sources down to K = 17.5 and, more importantly, because the background 
source counts do not deviate from the expected behavior until K > 17 (see § 4.3.2). Using the 
coordinates measured by Barsony et al. (1997, hereafter BKLT) for stars above their detection 
limit (K < 14.5) and below our saturation limit {K > 10), we derived relatively precise (±1") 
coordinates for the fainter sources {K > 14.5). 



3. Individual Source Characteristics 

3.1. Spectral Types 

The spectral types derived for the p Oph sample are listed in Table 1, along with measurements 
of Br7 emission and continuum veiling at 2.2 /xm (rx = /2.2(IR excess) //2. 2 (star)). Typical 
uncertainties are 0.25-0.5 in rx, depending on the signal-to-noise, spectral type, and degree of 
veiling. Most of the objects in our sample have been classified by the metallic absorption lines 
appearing in the X-band spectra, as described in § p.l.l| . In our X-band sample, we include all 
additional data obtained by Williams et al. (1995) (GY5, GYIO, GY37, GY202, 162618-242611). 
A total of 40 sources from this combined spectroscopic sample fall within the cloud core. We have 
adopted the spectral type of WGM for one of these (162618-242416) and now refer to it as part of 



the WGM sample. In § 3.1.2, we derive spectral types from the steam absorption in these spectra 
and compare our results to those of WGM. We revise the classifications of WGM for objects not 
in our sample to correct for a possible systematic error due to gravity dependence of the steam 
band strength. Comments on the individual sources are given in § The optical classification of 
the brown dwarf GY141 (Luhman et al. 1997) is also included in Table 1. 



3.1.1. IR Spectral Classification of the New Sample 

In a study of the small stellar aggregate in L1495E in Taurus, LR developed a technique of 
X-band spectral classification {R ~ 1000) to derive the spectral types and continuum veilings of 
young, late-type stars (~ 1 Myr, >G0). This method was then used to classify ~ 100 stars in 
IC 348, where the optical and IR spectral types were in reasonable agreement. The subtleties of 
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K-hand classification of young stars are discussed by LR and LRLL. For instance, the surface 
gravity of pre-main-sequence stars is intermediate between that of dwarfs and giants, resulting 
in weakened Na absorption at both 2.2 /xm and 8200 A in young M stars relative to M dwarfs 
(Luhman et al. 1998a; Martin, Rebolo, & Zapatero Osorio 1996). Comparing the data obtained in 
IC 348 {R = 800) and that presented here {R = 1200), we find that the higher resolution provides 
more accurate spectral types for G through M stars while maintaining sensitivity to faint sources. 
At higher resolution we obtain more accurate measurements of Mg at 2.28 fj,m, as illustrated in 
the observations of GY240A, GY240B, and GY250 at R = 800 and 1200 (Figs. |-D. The strength 
of Mg relative to Ca is important in the spectral classification of K through mid-M stars. Mg 
becomes weaker with later spectral type and eventually disappears beyond M2. For mid-M and 
later, two other features (Ca and an unidentified line) appear near 2.28 and can be mistaken 
for Mg when blended together at R < 1000. The two lines are resolved at R > 1200, allowing a 
more accurate classification (see V410 Tau X-ray 6 in Figure ^). 

The new spectra are shown in Figs. |2[|6| in order of spectral type and veiling. Background 
giants have very distinct spectra and are relatively easy to classify in moderate signal-to-noise 
data. A few spectra of lower quality are also shown, where the strong CO absorption could 
arise in either background giants or late-type cluster members. The spectral slopes, atomic line 
strengths, and mid-IR emission excesses imply that most of these objects are probably the latter. 
Given the various effects of veiling, surface gravity, and temperature on the appearance of these 
spectra, degeneracies in the classifications can occur. For instance, the sources classified as M4 
are also matched by late-K stars with atomic lines weakened by veiling and CO enhanced by low 
gravity. After comparison of IR and optical types in IC 348, we find that generally the simplest IR 
classification is the appropriate one, rather than a spectral type that requires simultaneous veiling 
and strengthening of CO. Previous optical and IR spectral types listed in Table 1 agree well with 
our measurements. Exceptions include SR20 and GY410, which were previously classified as GO: 
(Bouvier & Appenzeller 1992) and M6 (Greene & Meyer 1995), while we find very good fits with 
K5-M2 and K6-M0, respectively. Differences between our classifications and those of WGM are 
discussed in the next section. 



3.1.2. Spectral Classification with H2O Absorption 

In stars later than M5, the atomic absorption lines weaken and are only useful in classification 
of high signal-to-noise spectra. To derive spectral types for faint young low-mass candidates, we 
must instead rely on absorption in H2O (1.7-2.15 /im, ^ 2.3 ^m) and CO (> 2.29 /im). Because 
these features are broad and become stronger with cooler temperatures, lower spectral resolution 
{R < 500) is adequate for measuring spectral types. The CO absorption bands of young stars 
generally appear to arise from photospheres rather than disks (Casali &: Eiroa 1996; Greene & 
Lada 1997) and CO and steam should form in the same region of a disk, thus we expect the 
measured steam absorption to be dominated by stellar photospheres. 
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Detectable K-band steam absorption begins in early M dwarfs, causing the stellar continuum 
to peak at ~ 2.25 /Ltm and slope down to shorter and longer wavelengths. These depressions 
become more pronounced with later M dwarfs until M8V (WGM). Tokunaga Sz Kobayashi (1999) 
found that the i^-band spectra of the field dwarfs vB 10 (M8V) and DENIS-P J1058.7-1548 
(L3V, Kirkpatrick et al. 1999) are quite similar despite their very different optical spectra. It is 
likely that the steam absorption saturates for a large range of temperatures and is thus a precise 
indicator of spectral type only for <M8V or late L. 

Steam absorption can be measured through a reddening-independent index, as demonstrated 
by WGM. To estimate spectral types from such an index, the stellar continuum across the entire 
K-hand must be reliably measured, requiring observations of telluric standards immediately 
before and after the targets and the use of the same instrument configuration for all observations. 
Following these guidelines in observations with FSpec, we find that the spectral slopes are 
reproducible for the same object at different times, while there can be systematic differences from 
instrument to instrument (e.g., FSpec and Kitt Peak IR Cryogenic Spectrometer). 

For the classification of brown dwarf candidates in p Oph, WGM used a relation between 
steam absorption and spectral type observed for M dwarfs. However, since the surface gravities of 
these pre-main-sequence objects are ~ 2 dex lower than those of M dwarfs, we must investigate 
whether steam absorption is sensitive to surface gravity over the regime from dwarfs to young stars. 
Thus, we have compared the i^-band spectra of M dwarfs and optically classified young latc-type 
objects. Spectra covering 2.0-2.4 //m are provided in LR for the young objects V410 X-ray 3 
(M6, Strom & Strom 1994; Luhman et al. 1998a), V410 X-ray 5a (M5-M5.5, Strom &: Strom 
1994; Bricefio et al. 1998), and V410 X-ray 6 (M5, Strom & Strom 1994; M5.5, Luhman 1999). 
When comparing these data to a spectrum of GL 406 (M6V; same star used to classify Vx3 in 
the optical by Luhman et al. 1998a), we find that the steam in the young sources is stronger than 
in the dwarf counterparts. Therefore, whereas steam in dwarfs probably saturates near M8V, 
saturation may begin at M6 or M7 in young sources. The deep steam absorption present in young 
objects may seem contrary to the weak absorption observed at the very low surface gravities 
(log g = 0-1) of cool giants (Kleinmann &; Hall 1986). However, the synthetic spectra of Allard et 
al. (1995) do predict the strengthening of steam that we observe from M dwarfs (log g ~ 5.5) to 
pre-main-sequence objects (log g ^ 3.5). 

Using a steam absorption spectral sequence derived for evolved M dwarfs, the steam indices 
would imply spectral types of M6 for Vx5a and M7-M8 for Vx3 and Vx6, which are 1-2 subclasses 
later than their optical spectral types. To derive consistent spectral types from the steam 
absorption, the spectral sequence needs to be based on optically classified young stars. With Vx3, 
Vx5a, and Vx6 as spectral standards, we measured spectral types from the data of Williams et 
al. (1995) (GY5, GYIO, GY37, GY202, 162618-242611) and for the coolest objects in our new 
sample (GY29, GY59). As seen in Table 1, these spectral types are systematically earlier than 
those reported by WGM. Note that the M6-M7 spectral types may be upper limits because of the 
onset of saturation in steam. We classified an additional star in the WGM sample, GY84, with 
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the atomic absorption lines in the manner discussed in the previous section and used by LR and 
LRLL. We find a spectral type of early M for this object, in contrast to the M6 classification of 
WGM. 

With the later WGM spectral types, some of the sources are very cool yet luminous and fall 
high above the youngest theoretical isochrones. This placement appears to be inconsistent with 
photometric modeling (CRCTL). These sources are moved to the left on the H-R diagram with 
our classifications and have more plausible temperatures given their luminosities and likely ages. 
The revised parameters are more consistent with the photometry of CRCTL. 

The above analysis implies that the WGM spectral types are too cool by ~ 1 subclass for 
M6-M7 and ~ 1-2 subclasses for >M8. For WGM sources that were not in our sample or could not 
be classified in our data, we have revised the WGM spectral types accordingly, with the results 
of CRCTL as an added constraint. For 162618-242416, we adopt M5 as the spectral type, which 
is consistent with the data of both CRCTL and WGM. Given the apparent systematic error in 
WGM spectral types, the WGM classification of M8 for GY64 is revised to M7, which agrees with 
CRCTL. As noted by WGM, the JHK photometry of GYll indicates heavy continuum veiling, 
which would lead to an overestimate of the temperature when steam is used in the classification. 
Therefore, the WGM spectral type of M6.5 probably does not need to be adjusted to an earlier 
type, as suggested by the temperature estimated by CRCTL, which is equivalent to M7.5. Finally, 
three sources do not fall in our classified sample or the data of CRCTL. The relatively early 
WGM spectral types of GY31 (M5.5) and GY326 (M4) are adopted while the M8.5 spectral type 
of GY310 is revised to M7 for use in our analysis. 

3.1.3. A Possible FU Ori Object 

On the northeast outskirts of the cloud, the bright IR source known as IRAS64A (Wilking, 
Lada, &: Young 1989; Green et al. 1994) does not appear to be a typical pre-main-sequence 
star. The X-band spectrum exhibits weak atomic lines and strong CO, indicative of a K giant. 
However, the near- and mid-IR excess emission towards this source implies youth and cluster 
membership. As discussed in more detail in § we find that the data matches that of FU Ori 
objects (Hartmann & Kenyon 1996). However, this object may also be a background source, such 
as a post-AGB star or Mira variable. 

3.2. Extinction 

The recent survey by BKLT provides a homogeneous set of coordinates and JHK photometry 
for virtually the entire spectroscopic sample (see Table 1). The photometry is in the SQHD 
photometric system, which is similar to CIT (Kenyon, Lada, & Barsony 1998b). Photometry 
and coordinates for ROXs39 and SR22 are taken from Greene et al. (1994) and references 
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therein. The measurements of these two stars were converted to CIT using transformations from 
Humphreys, Jones, & Sitko (1984) and Leggett (1992). Since the p Oph cloud core represents 
a relatively small fraction of the square degree observed by BKLT, most of their survey region 
has low reddening and background stars are likely to dominate the stellar counts. Under this 
assumption, Kenyon et al. (1998b) used the IR color-color diagram to measure an extinction slope 
of E{J — H) / [H — K) = 1.57 ±0.03. We have spectroscopically identified several background stars, 
whose colors in Figure |7| are consistent with this reddening vector. For dereddening individual 
sources, we used the extinction law of Rieke & Lebofsky (1985). Although it differs from the slope 
found by BKLT, the difference is minor and the slope of BKLT alone cannot be interpreted in 
terms of individual color reddening coefficients. 

Figure ^ shows H — K versus J — H diagrams where the main sequence and the locus of 
classical T Tauri stars (CTTS) (Meyer et al. 1997) are plotted with their respective reddening 
bands. In the first three panels we present the spectroscopic sample as a function of r^- The 
sources behave as expected. Those with heavily veiled spectra and large H — K excesses fall on 
the reddened CTTS locus or to the right of it. Moderately veiled objects tend to fall between 
the reddened main sequence and the extension provided by the CTTS locus. Objects with no 
detectable veiling and normal colors are consistent with reddened main sequence stars. This 
behavior is similar to the trend observed by Casali & Matthews (1992) and Casali & Eiroa (1996) 
in spectra of CO absorption in a smaller sample of young stars. The scatter in this correlation, 
where a few sources with large rx have small excesses va. H — oi vice versa, can be explained 
by variability in the excess continuum emission (mostly at K). Such variability is evident in the 
data for GSS26 and GY51 (see § |A|.) 



Extinctions were estimated for stars earlier than GO by dereddening to the main sequence 
value of J — i/ for a given spectral type. Reddenings for late-type stars are often also estimated 
by assuming main sequence values for the intrinsic colors. In an older cluster such as IC 348 
(0.5-10 Myr), since much of the population consists of evolved Class III sources that show little 
evidence for disks, this approach for low-mass stars is reasonable. On the other hand, in a very 
young cluster such as p Oph (< 1 Myr), many of the stars are Class I through Class II with IR 
colors affected substantially by emission from circumstellar disks and envelopes. 

By using E{R — I) to deredden the J — H and H — K colors, Meyer et al. (1997) derived 
a locus of intrinsic colors for CTTS in Taurus, which are equivalent to Class II sources. The 
CTTS locus is shown as the dashed line in Figs. |^ and ^ and has been reproduced by models of 
star-disk systems (Lada & Adams 1992; Meyer, Calvet, & Hillenbrand 1997). The CTTS locus 
of Meyer et al. (1997) was measured for stars of spectral types near MO, where the origin of the 
locus falls at the main sequence colors of MO. As predicted by Meyer et al. (1997) and observed 
by Luhman (1999) in IC 348, the locus and its origin are shifted to redder H — K for later M 
spectral types. For earlier stars that fall within the reddening band of the CTTS locus, extinctions 
were calculated by dereddening the J — H and H — K colors to this locus. Sources to the right 
of the CTTS reddening band in Figure are heavily veiled and have uncertain spectral types. 
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consequently they are not plotted in the H-R diagram. 



3.3. Effective Temperatures and Luminosities 

To convert spectral types and photometry to effective temperatures and bolometric 
luminosities, we have adopted the temperature scale and procedures discussed in detail by LR 
and LRLL. We combined the dereddened J photometry with the bolometric corrections compiled 
by Kenyon & Hartmann (1995) and a distance modulus of 6.1 (Whittet 1974) to arrive at the 
bolometric luminosities. The source extinctions at J, effective temperatures, and bolometric 
luminosities are listed in Table 1, where the interstellar reddening law of Rieke & Lebofsky (1985) 
is used with Aj in the Johnson-Glass photometric system. 

No J photometry is available for two sources in our sample, GY244 and GY269. Since both 
show strong emission in the mid-IR, the H — K color may be contaminated significantly by excess 
at K. Our measurements of rx from the spectra could be used to correct for this effect, but since 
veiling can be variable and the spectra and colors were not measured simultaneously, we refrain 
from this approach. Although we cannot reliably place it on the H-R diagram, we will use the 
spectral type of GY244 and a canonical age to arrive at an approximate mass that can be added 



to the IMF for the cloud core in § 4.3.1. Such an estimate was not necessary for GY269, which 
falls outside of the cloud core. For the B star SI, we refer to the calorimetric luminosity derived 
by Wilking, Lada, & Young (1989). 



4. The p Oph Stellar Population 

4.1. Cluster Membership 

Classic methods of determining cluster membership (e.g., proper motions, Li absorption) 
are unavailable for most of the stars in our spectroscopic sample. However, we can confidently 
distinguish cluster members from field contaminants - foreground and background stars - through 
other means. With iiT-band spectra of modest signal-to-noise, background giants are easily 
identified through their distinctive continuum structure and strong CO absorption. All of the 
background giants in our sample appear in the less reddened area surrounding the p Oph cloud 
core, as indicated in Figure ||. The sample of WGM includes one likely background giant and two 
background early-type stars. GY297 shows little or no reddening in the IR colors and has no other 
indications of youth or cluster membership, such as near- or mid-IR excess emission or weakening 
of Na expected for young M stars (Luhman et al. 1998a, 1998b). It is a probable foreground 
star and is omitted from further analysis. WSB45 and WSB46 are a close pair (10") that show 
no obvious reddening in their IR colors or Li absorption (WSB45, Bouvier & Appenzeller 1992). 
Strong, variable X-ray emission and weak Ha emission (3-6 A) have been detected towards these 



- 11 - 



sources (Montmerle et al. 1983), suggesting they are either pre-main-sequence stars or active field 
dwarfs. We omit these stars for our discussion here. The remaining sources are unhkely to be 
background or foreground stars since they fall above the main sequence and have reddened colors. 
To examine cluster membership at fainter magnitudes beyond the limit of our spectroscopic 
sample, we will use a newly acquired KLF and diagrams oi H — K versus K to develop a model 
for the background star population shining through the cloud core (§ [4.3.5 ). 



4.2. The H-R Diagram and Star Formation History 



We use the estimates of T^s and Lboi from § 3.3 to place the il'-band spectroscopic sample on 



H-R diagrams in Figure All additional sources observed by WGM and Luhman et al. (1997) 
are also included. We have omitted objects with uncertain spectral types (<F8, K-M, featureless) 
and the foreground and background stars. The diagram of H — K versus K in Figure ^ indicates 
that the core spectroscopic sample is complete to -f^dereddened = llfori? — i^<2 {Ak < 3), 
corresponding to the dashed line on the H-R diagram for the core. 

To interpret the H-R diagram of p Oph, several sets of low-mass evolutionary tracks are 
available. Luhman (1998), LR, and LRLL made comparisons of the D'Antona & Mazzitelli (1994; 
DM94) and D'Antona & MazziteUi (1997; DM97) model predictions with observations of YY Gem 
and CM Dra, Pleiades, the main sequence, and globular clusters. Whereas calculations of both 
DM94 and DM97 were used in the analysis of data for IC 348 to provide continuity with previous 
studies, we will proceed with only DM97 in the following discussion. 

As shown in Figure ^, the tracks of DM97 imply a median age of ^ 0.3 Myr for the cloud 
core, with a few stars older than 1 Myr, particularly at the lowest masses where the scatter in 
apparent ages is largest. A similar range of ages is found in the sample surrounding the cloud core, 
where the median age is slightly older. A trend towards older ages with increasing distance from 
the cloud core is more obvious when these results are compared to the H-R diagram plotted by 
Greene & Meyer (1995) for optically visible stars observed across three square degrees (Bouvier &: 
Appenzeller 1992). Greene & Meyer also derived an H-R diagram for a sample observed with IR 
spectroscopy, composed of 5 stars from the core and 14 stars from the outer cluster. Using DM94, 
they found a distribution of ages very similar to that of the cloud core in Figure ^, although 
the study was slightly deficient in stars older than 1 Myr relative to the work presented here. 
Considering the uncertainties in model predictions and luminosity estimates at such early stages 
of evolution, absolute ages for these individual stars are not very meaningful. However, when the 
same set of evolutionary tracks is used from cluster to cluster, it is clear that p Oph is one of the 
youngest nearby stellar populations. 



4.3. The Initial Mass Function 
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4-3.1. The Spectroscopic Sample 

In conjunction with the evolutionary tracks of DM97, we can use Teg and Lboi to estimate 
masses for individual sources in p Oph and construct a cluster IMF. Since the cloud core provides 
an excellent screen against background stars, cluster membership and completeness are more 
readily addressed and we thus restrict the analysis of the mass function to sources within this 
region (see Figure ||). Our spectroscopic sample contains 39 objects that fall within the cloud 
core. Seven sources (GY6, GY20B, GY205, GY254, GY265, GY214, GY227) exhibit heavily 
veiled, featureless spectra and are excluded from this calculation of the IMF. Photometry for 
the companion GY20B is not available, hence it is absent from the figures. Three sources in the 
cloud core, GY192, 162712-243449, and GY192, had relatively uncertain late-type classifications 
and were not plotted in the H-R diagrams. They are indicated by circled dots in Figs. Il^ and 



12. GY192 has strong emission in IRAS and K — L' measurements and moderate excess in the 
J — H and H — K colors. We find that the spectral type (>M3) and luminosity estimate are 
consistent with the age of the cluster, implying a mass of ~ 0.2 Mq. Strong excess emission in 
H — K and K — L' and weak i^-band features are found in 162712 — 243449. The spectrum 
of GY128 is very noisy and quantitative measurements of its features are not possible. It also 
exhibits a flat-spectrum SED and large K — L' excess. Since 162712 — 243449 and GY128 cannot 
be classified and appear to be high-excess, heavily veiled objects similar to the seven featureless 
sources mentioned previously, they are omitted from the IMF. Two sources in our sample have 
uncertain reddenings and luminosities due to strong IR excess emission (GY51) and a lack of J 
photometry (GY244). By assuming these objects have ages comparable to the rest of the cloud 
core, we combined their spectral types with the evolutionary tracks to arrive at mass estimates. 
Because of uncertainties in the spectral types and subsequent mass estimates for two of the three 
early-type stars, the two highest mass bins are given widths of Alog M = 0.4. For GY182/WL16 
we use the -ff-band spectral type and a rough luminosity estimate to place it in the mass bin 
centered at log M = 0.35 (1.4-3.5 Mq) (Biscaya et al. 1999). For the remaining 23 stars from our 
sample and 6 stars from WGM and Luhman et al. (1997) shown in the H-R diagram of the cloud 
core, individual masses were estimated with the evolutionary tracks. The two lowest mass bins 
were also doubled in width to compensate for the greater uncertainties in the evolutionary tracks 
and temperature scales in this regime. 

The resulting IMF for the cloud core is composed of 36 sources, represented by the solid 



histogram in Figure 11. 



4-3.2. Completeness Correction for if < 13 

A completeness correction can be applied to the above IMF by 1) developing a model for 
the background stars shining through the cloud core at K, 2) identifying likely cluster members 
in the cloud core KLF that lack spectroscopy, and 3) estimating their masses by combining a 
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canonical cluster age implied by the H-R diagram with photometry and evolutionary tracks. We 
have applied this approach previously in our studies of L1495E and IC 348 (LR and LRLL). 

In the case of p Oph, we have obtained deep K-hand images of the cloud core {Ay ^ 50), 
covering the same region surveyed by CRBR (see Figure |l]). For the following discussion, we 
have merged our photometry with that of BKLT and CRBR. We use the data of BKLT when 
both H and K are provided {H ^ 15.5, K ^ 15.0), photometry of CRBR for the remaining 
sources detected in that work [H ^ 17, K ^ 15.5), and our new data for the faintest stars. The 
resulting set of photometry includes objects found in our survey and in that of BKLT but not 
reported by CRBR: 162706 - 243811, 162713 - 243330, 162726 - 244045, 162648 - 242836, and 
162622 — 242254 (GY12). The first four sources are faint companions to brighter objects while 
GY12 falls within bright nebulosity, explaining their absence from the photometry of CRBR. 
We reject five stars detected by BKLT that do not appear in our images or those of CRBR: 
162718 - 243433, 162720 - 243820, 162703 - 243726, 162719 - 244156, and 162719 - 244122. 

The cloud core KLF for the combined BKLT and CRBR data is given in Figure 0, m 
addition to the deeper extension provided by our data. It is evident from the rapid rise in counts 
at > 16 that the new KLF reaches the background star population behind portions of the 
cloud core. To simulate this field population, we selected all BKLT sources within a 25' x 10' 
area centered at a = 16''24™21!30, 6 = -24°42'00" (1950), which roughly equals the size of the 
cloud core. This position likely contains only relatively unreddened field stars, because it is near 
the outskirts of the BKLT survey where the molecular column density associated with p Oph is 
quite low. This off-held KLF must be reddened by Ak ~ 3.75 to fit the background star KLF 
(K > 16) observed towards the core. Such an average extinction is also consistent with the lack of 
background stars at i^' < 12 in the spectroscopic sample for the cloud core. The off-field KLF is 
fit well with a function of the form N[m)dm oc 10"™ dm, where a = 0.31, which is consistent with 
the Galactic field distribution observed by Wainscoat et al. (1992). When a KLF of this form is 
reddened, Comeron, Rieke, & Neuhauser (1999) found that the exponential form and value of a 
are preserved regardless of variable extinction or clumps. Thus, we show this fit to the background 



stars within the cloud core KLF in Figure 12. In a diagram of H — K versus K in the third panel 



of Figure |10|, we have reddened the off field by the best fit extinction derived in § [4.3.5| , illustrating 
where background stars should appear in the photometry. 

Nine stars showing either little or no features and high excess emission were excluded from 



the IMF determination in § [4.3. 1| since we could not constrain their spectral types or place them 
on the H-R diagram. Similarly, we do not attempt to add the high-excess, featureless objects to 
the IMF in the following completeness correction. Since many of the stars with no spectra at 
K < 13 have L'-band measurements, we will reject objects with K — L' ^ 2, which is typical of 
Class I sources. The mass function derived from this analysis will be based on the assumption 
that these objects are not biased towards one mass regime. This assumption is consistent with 
recent optical spectroscopy of Class I objects in Taurus, which exhibit a range of spectral types 
similar to that found among T Tauri stars (Kenyon et al. 1998a). 
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CRCTL detected mid-IR emission in ISO observations of low-mass candidates in the cloud 
core, several of which have been classified by WGM and Luhman et al. (1997) and added to the 



IMF in § 4.3.1. The four remaining sources that lack spectroscopy show weak mid-IR excess 
emission and can be added to the IMF by the photometric modeling described below. These 
objects (162623 - 242603, 162653 - 243236, 162710 - 242913, CRBR33) are represented by the 
open stars in the Figs. ^ and 12. 

The simulated background KLF predicts that only three stars in the K = 12-13 magnitude 
range are background. The cloud core KLF contains a number of sources in excess over the 
background model for K = 11.5-13. In addition, many of the objects between K = 12-13 have 
H — K > 3, much larger than expected for background stars, as shown in Figure |l^. We therefore 
assume that all the objects with K = 12-13 are members of the embedded population, except for 
three we rejected that have plausible colors to be candidate background objects. Within these 19 
remaining objects, 9 show evidence for strong IR excess emission in H — K oi K — L' and are 
excluded from the IMF modeling. 

The remaining ten sources are added to the IMF as a completeness correction in the following 



manner. As in § 3.2, extinctions are calculated by dereddening the J — H and H — K colors 
to the CTTS locus. When J photometry is not available, an intrinsic color oi H — K ^ 1 \s 
assumed, which is the maximum expected on the CTTS locus. Given the star formation history 
implied by the DM97 tracks in Figure P, we assume canonical ages of 0.5 and 1 Myr for these 10 
objects and convert the dereddened J (or H) magnitudes to masses for each age. After adding 
these photometrically derived masses to the IMF for the spectroscopic sample, we arrive at the 



distributions indicated by dashed lines in Figure 11. The completeness corrections for ages of 0.5 



and 1 Myr are very similar. With the exception of the sources with large IR excess emission and 
featureless spectra omitted from the IMF, the final mass function for the cloud core is complete to 
~ 0.08 Mq. 

The derived IMF is only the primary star mass function. Although a large number of binary 
systems have been resolved (separation > 0'.'005) in lunar occultation observations of p Oph 
(Simon et al. 1995), our derivations include only companions that are detected in the images 
of BKLT and CRBR (separation ^ 1"). Our IMF can be compared directly to results in other 
clusters, which generally also apply only to primary stars. As discussed by LRLL, under a broad 
range of assumptions, the inclusion of binary companions in the IMF produces a single star mass 
function with a low-mass slope of ~ 0.2-0.7 greater than that of the primary star mass function. 



4.3.3. General Behavior of the IMF 

At masses higher than ~ 0.4 M©, the DM97 IMF for p Oph matches that of Miller & Scalo 
(1979), while falhng less steeply than that of Scalo (1986). After the peak at ~ 0.4 M©, the IMF 
slowly declines to the hydrogen burning limit with a slope of ~ —0.5 in logarithmic units (where 
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Salpeter is +1.35). The lower limits for counts of young brown dwarfs suggest that the slope is 
not significantly less than —0.5 from 0.02 to 0.4 Mq. Allowing for the errors and incompleteness 
below the hydrogen burning limit, these results arc consistent with the logarithmically flat IMFs 
between 0.05 and 1 Mq estimated by CRBR and Strom et al. (1995) in IR luminosity function 
modeling of p Oph. 

Using the tracks of Baraffe et al. (1997), Bouvier et al. (1998) have recently reported a 
slope of —0.4 for the low-mass IMF in the Pleiades, which agrees with our results in p Oph and 
IC 348. LR, LRLL, and Scalo (1998) compared the IMFs measured near 0.1 Mq in photometric 
and spectroscopic studies of young, open, and globular clusters. The new results for p Oph and 
the Pleiades further support their suggestion that these various regions do not show substantial 
variation in the low-mass IMF. 

In comparison, recent discoveries of low-mass stars and brown dwarfs by 2MASS imply a 
present day field mass function that rises more rapidly with a slope of +0.3 (Rcid et al. 1999). 
In addition, there appear to be significant differences in the mass functions derived between 0.5 
and 2 Mq among these studies. The Pleiades IMF of Bouvier et al. (1998) is broad and peaks at 
0.4 Mq, consistent with our results in IC 348 and p Oph and the IMF of Miller & Scalo (1979). 
However, compared with this IMF, the population in Orion is deficient in stars over this mass 
range and is more sharply distributed at low masses, with a peak at 0.2 Mq (Hillenbrand 1997). 

Some of these variations may be due to differences in the techniques and tracks used in the 
young clusters. For instance, the Pleiades data were converted to masses with theoretical relations 
of M\^o[-mass (> 0.7 Mq, DM97) and Mj-mass (< 0.7 Mq, Baraffe et al.), whereas spectral types 
and H-R diagrams were used for Orion (DM94), IC 348 (DM94 and DM97), and p Oph (DM97). 
However, these explanations may not account for the behavior of Orion; whether it is different 
from other young clusters remains an open question. 

4-3.4- Comparison to the Pre-Stellar Clump Mass Function 

The proximity and extreme youth of the p Oph star forming region make possible the 
important comparison of the stellar IMF to the mass function of pre-stellar clumps. In deep 
millimeter continuum observations of p Oph, Motte, Andre, & Neri (1998) have studied emission 
associated with circumstcllar material, dense cores, and the ambient cloud. They detected 58 
starless clumps that showed no radio or IR sources and appeared to be pre-stellar or protostellar in 
nature. After estimating individual clump masses, they constructed a mass function that exhibited 
a slope of ~ 1.5 for ^ 1 Mq and ~ 0.5 for ^ 1 Mq. Motte et al. noted the similarity to the field 
star mass function (e.g.. Miller & Scalo 1979) and suggested a connection between the clump 
and stellar mass functions. Testi & Sargent (1998) have reached a similar conclusion in a study 
of Serpens. In this scenario, most of the mass in a given clump eventually accretes into a star 
or multiple system. With the new data presented here, we can now directly compare the stellar 
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IMF for the p Oph cloud core to the clump mass function measured by Motte et al. (1998), as 
illustrated in Figure |l^. Within the fairly low number statistics, the two mass functions are quite 
similar, providing further support for the notion that clump fragmentation directly influences the 
stellar IMF in a site of clustered star formation such as p Oph. 



4-3.5. Using Reddening Models to Test Membership at K > 13 
To explore the low-mass IMF further requires that we determine cluster membership for 



objects fainter than K = 13. Figure 10 shows that simple KLF modeling is inadequate for this 
goal. A number of faint {K > 13), relatively blue {H — K < 1.75) objects appear that may 
be low-luminosity cluster members rather than background stars. The isolation of these objects 
suggests that detailed modeling of the reddening that combines H — K colors with the KLF may 
be more successful in probing for very faint cluster members. 

There are discrepancies in the conclusions reached about the reddening in different studies. 
For example, CRBR confined their study to the region within the Ay > 50 contour in the CO 
maps of Wilking & Lada (1983). Based on a detailed discussion of the properties of the background 
reddened by this amount, they concluded that their sample included a number of background 
stars but not in a quantity that would significantly affect their results. In contrast, Kenyon et 
al. (1998b) used the H — K versus J — H color-color diagram to argue for an extinction only 
half as great as found from the CO maps, and they contend that CRBR may have sufficiently 
underestimated the contamination by background sources to affect their conclusions. Neither of 
these studies makes a quantitative determination of the effects of extinction variability on their 
conclusions. Variations within the CO beam would tend to increase the number of background 
sources seen through the cloud above the estimate by CRBR (as they point out). Variations 
will produce an underestimate of the average extinction in the color-color plot of Kenyon et al. 
(1998b), unless great care is taken to avoid selection effects against the redder sources. 

With the availability of spectra and ISO data for many of the faint objects in the CRBR 
study (CRCTL; WGM; this paper), it appears indeed that most are bona fide embedded 
cluster members, in agreement with CRBR rather than Kenyon et al. (1998b). We still need to 
reconcile the derivation of strong background contamination by Kenyon et al. (1998b) and their 
determination of an extinction level less than half of that from CO maps and a "normal" ratio of 
CO column and Ay (Wilking & Lada 1983). One contributing factor is that BKLT and Kenyon 
et al. (1998b) included more area outside the Ay = 50 contour in their study than did CRBR. 
Another possibility is that extinction variations may have significantly affected the near-IR color 
distributions. 

To pursue the latter question, we have used a simple model assuming a gaussian distribution of 
extinctions. We have confined our study to the same area studied by CRBR. Most of this area lies 
within the CO-derived Ay = 50 contour. Our model can be used to take a background population 
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and redden it in a manner that agrees with the observed properties of the background seen to 
be penetrating the p Oph obscuring cloud in our deep H and K images. From the background 
population appearing in off-cloud fields of BKLT, the cloud core KLF, and an assumption of 
uniform extinction, we had deduced Ak ~ 3.75 through the cloud (§ |4.3.2 ). However, in deep 



//-band images of selected regions within the cloud core, we found that the background begins to 
emerge aX H — K ^ 1.7b. This color is bluer than expected for the average level of Ak ~ 3.75, 
indicating that the extinction is non-uniform, so surveys at progressively shorter wavelengths are 
biased toward progressively less obscured members of the background. 

We have two independent measurements to constrain the behavior of the extinction. First, 
by comparing the on-cloud and off-cloud star counts at K, we have determined the net effective 
extinction at that wavelength. Second, by imaging deep enough at both H and K to see the 
background population through the cloud, we have a measure of the bias discussed above. The 
interpretation of this bias is relatively simple at these colors, since plausible background stars all 
have nearly the same H — K colors. We have assumed a uniform value of 0.15, which should be 
within about ±0.1 of the color of most background stars (e.g., Leggett 1992) 

We have modeled the extinction of the background assuming that the variations are Gaussian 
around an average value. We have taken a net Ak = 3.75 as a lower limit, to allow for the 
possibility that the "off' field we have used for comparison with the on-cloud field has some low 
level of extinction. We have assumed an upper limit of the net Ak < 4.05. From comparing the 
H and K luminosity functions, we find that the color difference at which the background emerges 
is 1.7 < H - K < 2.2. 

A range of parameters can fit the observed behavior within these limits. The true average 
extinction could lie from Ak = 3.9 with a standard deviation of 0.6 to Ak = 5.0 with a standard 
deviation of 1.45. A representative fit is Ak = 4.4 (i.e.. Ay ~ 40) with a standard deviation 
of 1.2 (see Figure p^ ), which is our best estimate of the true behavior of the cloud. Modeling 
the extinction variations has resulted in an increase by 10 to 20% in the estimated average level. 
However, the possibility that the CO measures are slightly biased toward a high estimate is also 
confirmed. 

Our extinction model also allows us to estimate the numbers of background sources we might 
see within our sampled region, as a function oi H — K (and K magnitude). The proportion is 
almost independent of K magnitude. The models with the largest background contamination 
are those with high average extinction, since the two-parameter fits then require large standard 
deviations for the extinction. For these "worst case" models, ~ 16% of the background stars would 
have H — K < 1.5, roughly independent of magnitude, and ~ 9% would have H — K < 1.25. Most 
of the successful models predicted signficantly lower levels of blue background star contamination; 
the values quoted are upper limits. These predictions are used in the following section to examine 
the status (reddened background stars versus low-luminosity cluster members) of the very faint 
{K > 14), "blue" {H-K < 1.5) sources. 
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4-3.6. New Suhstellar Candidates 

Figure [l^ shows that we cannot achieve complete identifications of cloud members from the 
relatively blue {H — K < 1.5) colors discussed above; at the brighter levels, many bona fide 
members are redder than this limit. However, we can use the predictions to see if we are sampling 
a portion of the very low-mass cloud members. For example, between K of 13 and 16, there are 25 
objects of which 9 are bluer than H — K = 1.5. Our background model would have predicted only 
4 objects that blue. The range if = 13 to 16 corresponds to masses from about 0.01 to 0.08 Mq, 
so the presence of a significant number of objects in this range is consistent with our previous 
arguments that the IMF includes a reasonably large number of brown dwarfs at least down to 
0.02 Mq. 

Because of the small numbers of sources included in our deep survey in these magnitude 
ranges, this conclusion is not well-established. Our main goal is to illustrate a method to extend 
cloud member identification to very low masses. One needs to image the cloud in at least two 
bands to a sufficient depth that the background stellar population is well characterized in both. 
The information on the background star color distribution allows construction of a model of the 
variations in extinction through the cloud. The modeling can be quite simple at H, K, and L 
because nearly all stars have similar intrinsic colors over these bands. Such a model can test 
whether relatively blue sources are likely to be cloud members as opposed to background sources 
shining through relatively thin sight lines of the cloud. 



4.4. SED Classes and the First Stages of Stellar Evolution 

Through low-resolution {R = 500) JHK spectroscopy, Greene k. Lada (1996) examined the 
near-IR properties of a large number of young stars in p Oph and other star forming regions. A 
range of IR SED classes was represented in their sample, including Class I (a > 0.3), flat spectrum 
(0.3 > a> —0.3), II (—0.3 > a > —1.6), and III (a < —1.6), where a is the mid-IR spectral index 
(Lada 1987; Greene et al. 1994). This index measures the slope of an SED from 2 /im to 10-20 /im 
and is sensitive to the presence of cool circumstellar material. Greene & Lada observed a trend of 
increasing line emission and reddening in addition to weakening absorption lines from Class III to 
Class I sources, consistent with the evolutionary interpretation of these observational classes. Our 
results on the relationship between veiling and location on the H — K versus J — H diagram shown 
in Figure ^ confirm these arguments. In addition, we find a decreasing incidence of detectable Br7 
emission with decreasing IR excess, as expected. 

By combining previous near- and mid-IR photometry with our measured spectral types, 
continuum veilings, and emission line strengths, we now investigate in more detail fundamental 
questions regarding the evolution of young stars. For instance, are some objects previously 
identified as Class I in fact older and partially emerged from their natal dust shells, but so heavily 
reddened by the molecular cloud that their SEDs place them in Class I? With these data, we will 
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1) deredden the spectral indices and review the SED classifications, 2) calculate the timescale for 

the Class I protostellar stage of evolution, 3) use our spectroscopic data to constrain the origin 
of IR excess emission in Class I and Class II objects, 4) estimate the dctcctability of absorption 
features in spectra of Class I sources, 5) discuss the nature of the reddening in Class II sources, 
and 6) compare the distribution of SED classes inside and outside of the p Oph cloud core. 

4.4- 1- Revised SED Classifications 

To sort the p Oph sources into SED classes requires a determination of their intrinsic 
K{2.2 /xm) — iV(10.6 pm) spectral indices. We have done so for a sample of all sources with 
K < 12 that have spectral types of Kl or later or exhibit featureless K-band spectra. The multiple 
systems GY20 and GY240 are not included since the components arc unresolved in the mid-IR 
data. IRAS64A, the possible background and foreground stars, and several sources with uncertain 
late-type classifications are also omitted. 

Ground-based mid-IR measurements were taken from Elias (1978), Lada & Wilking (1984), 
Young, Lada, & Wilking (1986), Wilking, Lada, & Young (1989), and Greene et al. (1994). Two 
objects were classified with the ISO observations of CRCTL. We avoided IRAS measurements 
because the beam is large enough that confusion with a fainter and redder object is a possibility. We 
have used these data with the values of Aj and K in Table 1 to compute the extinction-corrected 
K — N (2.2 — 10.6 ^m) spectral indices, a(der), given in Table 2. When available, a(2.2 — 20 /im) 
and measurements of continuum veiling and Br7 emission from the iC-band spectra are also listed. 

Stars within the reddening bands of the main sequence and CTTS locus should have 
reasonably accurate extinction estimates, so for them wc have relied on the dereddened indices in 
classifying their SEDs. When mid-IR data were available but extinctions could not be computed 
due to the lack of J-band photometry, the sources were classified by considering both a(obs) and 
tk- For sources to the right of the CTTS locus, scattered light may dominate and the reddenings 
are uncertain, thus we vaake. no attempt at dereddening their spectral indices. GY6, GYlll, 
GY214, GY254, GY265, GY269, GY274, and GY378 fah in this category and are likely true 
Class I rather than heavily reddened Class II objects given their strong veiling and steep SEDs 
between 2 and 20 pm.. GY205, GY224, and GY227 are flat-spectrum sources by their observed 
indices, which cannot be dereddened reliably due to large H — K excess (GY205) and lack of 
J photometry (GY224, GY227). The spectra of these three objects are featureless, supporting 
their flat-spectrum status, although the low signal-to-noise provides weak constraints on the 
veiling {rx > 0.5). Because GSS26, GY21, GY51, GY279, and GY315 are within the CTTS 
reddening band or a small distance to the right of it, we can calculate dereddened spectral indices. 
Although the dereddened indices for GSS26, GY21, and GY51 indicate Class II, they remain 
close to the fiat-spectrum regime, so we cannot confidently classify these as either fiat or Class II. 
Similarly, GY279 and GY315 are likely Class II, but we cannot rule out flat-spectrum status due 
to variability in GY279 and uncertain extinction in GY315. No mid-IR photometry is available for 
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GY81, but comparing the spectroscopic properties to those of the sources we have discussed, it is 
probably flat-spectrum or Class II. Other sources without mid-IR data are tentatively classified in 
a similar fashion. We have no J photometry for GY244, but the lack of veiling or Br7 emission 
implies little excess at K, hence we estimate the extinction from H — K. The dereddened spectral 
slope indicates that GY244 may be a heavily reddened Class II object. 

For calculating the distribution of sources across SED classes, we divide the 6 objects labeled 
as fiat/II equally between fiat-spectrum and Class II and consider flat-spectrum a part of Class I. 
With the revised classifications, there are 12.5, 26.5, and 33 sources among Class I, Class II, and 
Class III, respectively. The dereddening has moved 18-30% (8-13/44) of sources that are Class I or 
Class II by the observed indices from one class to another, in contrast to the suggestion of Greene 
et al. (1994) that less than 10% of sources should change SED classes upon dereddening. These 
discrepant conclusions may be caused by differences in how the dereddened indices were calculated 
in the two studies. Greene et al. state that a(2-10 jim) changes by 0.5 when reddened by Ay = 50 
under the extinction law of Rieke & Lebofsky (1985), whereas we find the change in the index is 
1.7. 

Selection effects will influence these results. For example, an object of a given mass is 
expected to fade in bolometric luminosity between Class I and Class III, favoring the detection of 
the Class I stage. However, it is also possible that a smaller portion of the luminosity emerges 
in the near-IR due to the larger extinction in the Class I phase, favoring the detection of the 
later classes in our K-ha,nA limited sample. For instance, 5/7 sources classified as Class I or 
flat-spectrum by Greene et al. (1994) have K > 12 and therefore are not included in the statistics 
we have discussed. Disks around very young objects where the dust has not started to coagulate 
or dissipate may hide the object from near-IR surveys along some sight lines. For example, if the 
mass of the planets in the solar system were distributed uniformly in a disk of 100 AU radius and 
3 AU thickness, the density would be of order 10~^^ g cm~^, corresponding to Aj ~ 1000 taking 
typical interstellar parameters. Unresolved binaries will tend to favor Class II designation. In the 
case of a Class II source paired with a Class III one, the mid-IR excess would bias the composite 
to Class II. A Class II or III source paired with a Class I one would produce composite near-IR 
characteristics (photosphcric absorptions, bluer K — N color) indicative of Class II. Nonetheless, 
the ratios of the numbers of sources in the various classes should give a rough idea of the true 
distribution among the SED classes. 

4-4- ^- Timescale for the Protostellar Stage 

A variety of arguments regarding the very early development of a young star assign probable 
durations to the stages corresponding to differing SED characteristics. Specifically, it is thought 
that the Class I SED type should persist for only 1 - 2 x 10^ years (e.g., Cohen & Kuhi 1979; 
Stabler 1983; Shu, Adams, & Lizano 1987). However, the photometric surveys of CRBR and 
Greene et al. (1994) have found Class I objects in larger numbers than this short duration would 
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suggest. 

However, after dereddening the sources and extracting statistics from a magnitude-limited 
sample, only about 17% (12.5/72) of them are Class I. This proportion is substantially lower than 
previous estimates (CRBR; Greene et al. 1994). If we interpret the proportion of Class I sources 
as an indication of the lifetime of that phase of stellar evolution and assume that star formation 
has occurred at a uniform rate over the last 0.5-1 Myr in the p Oph region, then the lifetime of 
the embedded stage is ^ 0.075 — 0.15 Myr. This estimate removes the discrepancy with theoretical 
expectations found in previous estimates of the proportion of Class I sources. 

4.4-3. Nature of Infrared Excess Emission in Classes I and II 

In a generally accepted view of the protostellar stage, the Class I SED has been modeled as a 
star surrounded by an accretion disk and a thick circumstellar envelope, where the disk is thought 
to contribute a substantial portion of the excess emission at 10 //m (Myers et al. 1987). Because of 
the lack of disk absorption features in Class I and II systems, Calvet, Hartmann, & Strom (1997) 
have suggested that the stellar magnetosphere induces a hole in the inner region of the disk, thus 
preventing the formation of a disk photosphere. They propose that a warm infalling envelope 
dominates circumstellar excess at 2 /xm. In a Class II object where the envelope is absent and 
only the disk remains, relatively little if-band continuum emission {rx ^ 1) is predicted (Meyer 
et al. 1997). We now consider this model in light of the new data we present. 

In May of 1996, GSS26 and GY51 (fiat/II) exhibited strong Br7 emission (2.9±0.2, 5.0±0.5 A) 
and veiling [vk ~ 4, 2), while in July of 1994 the Br7 emission was weaker (1.7±0.3, 1.5±0.5 A) 
and the excess was lower (r^ ^ 0.75 and 0.5). Furthermore, when these spectra are flux calibrated, 
we find that GSS26 and GY51 were ~ 1.2 and 1.3 mag brighter at K in 1996 relative to 1994, 
consistent with the increase in excess continuum emission implied by rx. Both the excess emission 
and Br7 line flux in each star brightened by factors of 5-10. If the Br7 emission arises in columns 
of material accreting from the disk onto the star (Muzerolle, Calvet, &; Hartmann 1998), then 
the simultaneous increase in the Br7 emission and the continuum veiling suggests that the excess 
emission is produced in a region directly associated with the accreting material (i.e., the disk) 
rather than a circumstellar envelope. 

We can also use the behavior of r^ with changing SED class to probe the origin of the excess 
emission. If an infalling envelope is responsible for the significant continuum veiling at 2 //m seen 
in a given source, strong mid- and far-IR emission should also be generated, appearing as a Class I 
SED. However, we find strong veiling, tk = 1-4, in several sources lacking this evidence of an 
envelope, as shown in Table 2. Examples include GSS26, GY21, GY51, GY279, and GY315, which 
are possibly Class II after dereddening. Even without dereddening, EL24 (rx = 1-3), GY129 
{rK > 1), GY167 (rx = 1-3), GY168 (r^ = 2), GY292 (rx = 1), and SR4 (rx = 1.5-2) are 
Class II, yet they have moderate to heavy veiling. Finally, we measure ri^: = 1 in GY262 and SR20, 
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which show no sign of an envelope by their Class III status. On the other hand, a disk is likely 
to be present in all of these SED classes. If the disk is the source of the K-hand excess emission, 
then the inner hole must be within the region where dust would emit strongly at 2 ^m while large 
enough to prevent dust destruction and the resulting formation of a substantial disk photosphere. 
The central hole in the model by Greene & Lada (1996) was tuned in this manner and produced 
veilings of rx ^ 1 for a pure re-processing disk and tk ^ 5 for a luminous disk, although this 
choice of hole size has been subsequently criticized by Calvet et al. (1997) as nonphysical. 



4-4-4- Detecting Absorption Features in Class I Sources 

In addition to Class II and Class III sources, can photospheric features also be detected 
in extremely red objects designated as Class I? Greene &: Lada (1997) observed absorption 
features in two sources, GY21 and GY279, and demonstrated that the lines likely originate in 
the photospheres of the central stars rather than disks. Although these objects were previously 
classified as flat-spectrum, they appear to be Class II when their spectral indices are dereddened 



(see § 4.4.1). With our large set of moderate-resolution spectra, we can further examine the 
detectability of absorption features in flat-spectrum and Class I sources. Featureless -fC-band 
spectra are found for GY6, GY20B, GY214, GY254, GY265, GY274, GY378 (Class I), GY205, 
GY224, GY227 (flat spectrum), and GY129 (Class II). However, absorption features are detected 
in the Class I objects GYlll and GY269, comprising about ~ 25% of Class I sources. Since the 
signal-to-noise is rather low in some of the featureless spectra, higher quality data could reveal 
absorption features in still more of the Class I stars. 



4.4-5. Nature of the Reddening in Class II Sources 

To test the origin of the obscuration for the highly reddened Class II objects, we have 
computed the average Aj for the Class II and Class III sources in the sample described here. The 
ratio (11)/^ J (III) = 1.1 lb 0.2; thus there is no excess of extinction for Class II objects. The 
distribution of extinction values is also similar between the classes. Thus, the Class II and III 
sources appear to be distributed throughout the cloud in a similar fashion, as one would expect. 
This result supports our conclusion that some sources previously considered to be of Class I or II 
are instead heavily reddened Class II or III. In addition, any circumstellar disks around Class II 
objects appear to have little influence on the extinction - either the disks are optically thin at J, 
or they are thick enough that systems oriented with their disks edge on are so heavily obscured 
that we do not detect them. 
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4-4-6- Spatial Distribution of SED Classes 

We now make a comparison of the IR characteristics of the areas separated by the dashed 
hne in Figure ^, defining stehar populations from the cloud core and the surrounding region. In 
the diagrams of H — K versus J — H in Figure ^, we find that a larger fraction of sources within 
the cloud core have near-IR excess emission. To be more quantitative, we designate sources with 
evidence for actively accreting circumstellar disks as those exhibiting Wx(Bt'j) > 1 A or r^^ > 0.5. 
In the core of p Oph, if we omit the B star and seven sources with uncertain veiling measurements, 
the frequency of disks is 17/27 (63%). The frequency is somewhat lower outside of the cloud core, 
23/45 (51%), after omitting the four earliest stars, two K-M stars, IRAS64A, and the foreground 
and background stars. The difference between the two populations is more pronounced when we 
examine the SEDs. For sources in the cloud core that have been classified in Table 2, the number 
distribution among Class I/flat, Class II, and Class III is 6.5, 9.5, and 10, respectively. Outside of 
the cloud core, on the other hand, the SEDs are generally more evolved, where the distribution is 
6, 17, and 23. 

In a study of the older cluster IC 348 (0.5-10 Myr) by LRLL, ~ 24% of sources within the 
core with ages less than 3 Myr met the W\(BT'y) > 1 A or > 0.5 criterion for active accretion. 
The lower fraction of disk sources in the core of IC 348 relative to the remainder of cluster was 
attributed to the accelerated truncation of disks at higher stellar densities (Herbig 1998; LRLL). 
In p Oph, since the densities are rather low throughout the region, the higher disk frequencies 
from the core outward are instead likely due to the birth of stars in the cloud core and their 
dispersal outward. Such an age gradient is consistent with the ages implied by the H-R diagram 



of the cloud core and surrounding area, as discussed in § 4.2 



5. Conclusion 



We have performed K-hand spectroscopy towards ~ 100 stars within the p Oph star forming 
region. We have measured spectral types and continuum veilings and have combined this 
information with previous mid-IR measurements and new near-IR imaging of the cloud core. Our 
conclusions are as follows: 



1. The K-h&nd steam absorption, known to strengthen rapidly with from mid-M to M8V and 
possibly saturate for later types, is also sensitive to surface gravity and differs between 

~ 1 Myr and the main sequence for a given spectral type. Thus, young late-type objects 
classified in the optical should be used as spectroscopic standards when using steam to 
measure spectral types for young, cool objects. 

2. Using the evolutionary tracks of DM97 to interpret the H-R diagram for p Oph, we estimate 
stellar ages ranging between 0.1 and 1 Myr, with a few stars that may be slightly older. 



-24- 



3. Prom a deep KLF of the cluster core we develop a reddening model for the background star 
population and identify likely cluster members falling below the limit of the spectroscopic 
sample. We use the age from the spectroscopic sample to convert their luminosities to 
masses. After the addition of these sources, the resulting IMF is complete to ~ 0.08 Mq. 
With the tracks of DM97, the IMF for p Oph matches that of Miller & Scalo (1979) at 
masses higher than 0.4 Mq. The IMF peaks at this mass and slowly declines to the hydrogen 
burning limit with a slope (in logarithmic units) of ~ —0.5, as compared to slopes of +1.35, 
0, and -2.6 for Salpeter (1955), Miller & Scalo (1979), and Scalo (1986), respectively. The 
exact shape of the mass function remains dependent on the theoretical evolutionary tracks 
and temperature scales. 

4. Our lower limits on the numbers of substellar objects demonstrate that the IMF probably 
continues a slow decline and does not fall precipitously below the hydrogen burning limit, at 
least down to ~ 0.02 Mq. 

5. The derived IMF is consistent with previous findings that the p Oph IMF is roughly flat 
from 0.05 to 1 M©. 

6. The stellar IMF in p Oph is qualitatively similar to the mass function of pre-stellar clumps, 
indicating that cloud fragmentation may play a direct role in determining the IMF. 

7. Very deep imaging at H and K (also at L) can be used to construct models of the extinction 
through the cloud and determine background contamination as a function of color and 
magnitude. We have used this method in a limited area to identify several faint sources 
{K = 13-16) that may be cluster members by their low reddenings {H — K < 2), and thus 
are extremely low-mass (10 to 80 Mj) substellar candidates. 

8. Extinction can significantly alter the true distribution of SED classes. The classifications 
of ~ 25% of sources originally classified as Class I-II are changed upon dcrcddcning of the 
mid-IR indices {a(2-10pm.)). After revising the SED classification accordingly, ~ 17% of the 
p Oph stars are Class I, implying ~ 0.075 — 0.15 Myr for the lifetime of this embedded stage. 

9. We detect absorption features in the spectra of ~ 25% of Class I and fiat-spectrum sources, 
supporting the feasibility of studying the photospheres of extremely young protostars. 

10. In observations of 20 sources separated by two years, we find that the continuum veiling 
and Br7 emission change simultaneously in two stars. Because the Br7 emission likely 
traces material in the accretion columns between the disk and the star, this behavior 
suggests the disk as the origin of the 2 pm. excess rather than an infalling envelope. The 
presence of significant veiling {rx = 1-4) in the spectra of several stars from flat spectrum to 
Class III also indicate that a disk, or at least something other than an envelope, is capable 
of producing large amounts of continuous 2 pm emission. 
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11. The IR properties of the sample imply a slightly more evolved population distributed around 
the cloud core, possibly due to the birth of stars within the core and their subsequent 
dispersal, which is also consistent with the ages implied by the H-R diagrams of the cloud 
core and surrounding region. 

We thank M. Rieke for helping to obtain the IR spectra and M. Meyer and I. N. Reid for 
providing results prior to publication. Detailed comments on the manuscript by C. Lada and M. 
Meyer were greatly appreciated. We are grateful to F. Allard, I. Baraffe, A. Burrows, and F. 
D'Antona for providing their most recent calculations and useful advice. This work was supported 
by NASA grant NAGW-4083 under the Origins of Solar Systems program. 

A. Notes on Individual Sources 

With the exception of RNOIC, which is similar to the data for the K2III k Oph (Kleinmann & 
Hall 1986), the ii'-band spectra of FU Ori objects are all similar and show characteristics distinct 
from normal stars, as shown in Figure |l^. Some of this difference is likely due to broadening 
of lines from disk rotation (Greene & Lada 1997). The FU Ori spectra also show two features 
on either side of where Mg appears in the giant (possibly the two lines seen in late M dwarfs) 
and unidentified absorption at 2.33 /xm. All of these characteristics are found in the spectrum of 
IRAS64A. While the FU Ori objects have strong CO absorption and late-type IR spectra, they 
tend to appear as F or G giants in the optical. This discrepancy is reconciled by attributing the 
CO to absorption in low-gravity disk atmospheres. However, we find that the optical spectrum 
(6000-9000 A) of IRAS64A is identical to the M8-M9III spectra presented by Kirkpatrick, Henry, 
&: Irwin (1997). High-resolution spectroscopy on the CO band heads is necessary to determine 
whether a circumstellar disk is the origin of these peculiar optical and IR spectra. 

Another unique source, GY182/WL16, is well-known for its strong CO band head emission, 
modeled successfully in terms of an inversion layer in a circumstellar disk (Carr et al. 1993; Najita 
et al. 1996). As seen in Figure ^, no photospheric features appear in the if -band spectrum of this 
Class I object. However, in observations of second overtone CO emission, Biscaya et al. (1999) 
detect absorption in the H-hand Brackett series of hydrogen, presumably arising from the stellar 
photosphere of the central star. They conclude that the spectral type of WL16 is between B8 
and A7 (Tefj = 8000-12000 K). Subsequent modeling of the CO emission further constrains the 
temperature of the star to be 9000 to 12000 K. 

GY20/DoAr24E is a binary system with a separation of 2'.'06 and component magnitudes of 
K = 6.9 and 7.9 (Simon et al. 1995). The curvature in the slope of the Class I secondary is due 
to its proximity to the primary and the resulting difficulties in extracting the spectrum. Since the 
companion is much redder than the primary, we assume the BKLT J photometry applies to only 
GY20A when deriving the luminosity. 
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For GY29, both the CO bands and the shape of the continuum produced by steam absorption 
match well with the spectrum of V410 X-ray 6 (M5.5) artificially reddened by Ay = 25. 

In addition to absorption in Br7, we clearly detect weak features of Na, Ca, and CO in 
the B star GY70/S1. The source of these metal lines is likely the fainter, unresolved secondary 
{AK = 1.8, 0'.'02, Simon et al. 1995). 

GY240/WL20 is a triple system (BKLT, references therein) appearing as one source in the 
data of BKLT and Wilking &: Lada (1989) and two sources separated by 2-3" in the images of 
CRBR and Barsony et al. (1989). We refer to the west and east components as GY240A and 
B. Since GY240B has no J photometry and the H — K colors of the two sources are similar, 
we derive the extinction from BKLT J — H for the composite system and assume it applies to 
each component. To estimate a J magnitude for the individual sources, we assume a difference of 
^ 1.2 mag, which is measured by CRBR at H and K. 

GY250/SR12 is a binary system with a separation of C/.'S and component magnitudes of 
K = 9.3 and 9.4 (Simon et al. 1995). We assume the components have the same spectral types 
and divide the luminosity in half. Only one source is added to the IMF, for reasons discussed in 

GY372/VSSG14 is an unresolved binary system (C'l, Simon et al. 1995) where the secondary 
is 0.8 mag fainter at K. As with GY70/S1, Br7 is strong while we also detect absorption in metal 
lines and CO that may arise in the secondary. Due to possible line and continuum veiling (from 
either the secondary or IR excess emission), the classification of the primary is uncertain (<F8). 
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Fig. 1. - All sources in p Oph observed with K-hand spectroscopy, with the exception of seven 
objects falling outside of this region (IRS2, IRAS64A, SR22, WSB38, WSB45, WSB46, ROXs39). 
Filled circles represent stars classified as background giants in this work and open circles are all 
other sources in our sample. Open triangles and the open square are the additional low-mass objects 
observed spectroscopically by WGM and Luhman et al. (1997). Filled triangles are identified as 
likely background stars by WGM. The area within the dashed line was surveyed through IR images 
by CRBR and corresponds to the cloud core {Ay ^ 50; Wilking & Lada 1983). Deeper iiT-band 
photometry of the same region was obtained for this work, in addition to J and H images towards 
the positions outhned by dotted lines. The origin corresponds to a = 16''24™07!7, S = -24°24'00" 
(1950). 



-31 - 



3.25 



2.75 



2.5 



2.25 



cd 
-i-j 

m 

o 
o 



fin 



1.75 



1.5 



1.25 



0.75 



— I — 1 — 1 — 1 — 1 — I — 1 — I 
Fe Mg Ti CO 



e NaTi 




2.15 



2.2 



2.25 



2.3 



2.35 2.15 



2.2 



2.25 



2.3 



2.35 



A (Aim) 



Fig. 2. — K-band spectra at i? = 1200 of early-type stars and sources exhibiting no detectable 
continuum veiling [tk < 0.25), with additional data at i? = 800 (*). Spectra are normalized at 
2.2 nia with constant offsets. 
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Fig. 3. — if-band spectra at R = 1200 of late-type sources exhibiting no detectable continuum 
veiling {tk < 0.25), with additional data at i? = 800 (*). Spectra are normalized at 2.2 /xm with 
constant offsets. Possible foreground stars (F) are also indicated. The IR spectrum of V410 Tau X- 
ray 6 (LR), which is classified in the optical as M5.5 (Luhman 1999), is shown for comparison. The 
spectra classified as M4 are also fit marginally by veiled late-K stars with enhanced CO absorption. 
Spectra are normalized at 2.2 /xm with constant offsets. 
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Fig. 4. — ii'-band spectra at R = 1200 of moderately veiled sources (0.25 < tk < 1), with 
additional data sX R = 800 (*). Spectra are normalized at 2.2 jim. with constant offsets. 
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Fig. 5. — if -band spectra at R = 1200 of heavily veiled sources {vk > 1) with additional data 
at i? = 800 (*). The curvature in the spectrum of GY20B is not real and is due to difficulties in 
extracting the spectrum separately from the companion GY20A. Spectra are normalized at 2.2 /xm 
with constant offsets. 
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Fig. 6. — i^-band spectra at R = 1200 of heavily veiled sources {vk > 1), background giants (B), 
and objects with uncertain late-type classifications, with additional data at R = 800 (*). Spectra 
are normalized at 2.2 fim with constant offsets. 




Fig. 7. — H — K vs. J — if for the ivT-band spectroscopic sample is given in a), b), and c) as 
a function of the continuum veihng, rx, measured in the spectra. Baclcground giants, early- type 
stars, and sources with uncertain late-type classifications are given in d). The main sequence 
(dotted line) and CTTS locus (dashed line) are shown along with the reddening bands of each, 
where the reddening vectors were derived from BKLT data by Kenyon et al. (1998b). 
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Fig. 8. — H — K vs. J — H for the i^-band spectroscopic sample within and outside the cloud core, 
as defined by dashed boundary in Figure ||. Background and foreground stars have been omitted. 
The main sequence (dotted line) and CTTS locus (dashed line) are shown along with the reddening 
bands of each, where the reddening vectors were derived from BKLT data by Kenyon et al. (1998b). 
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Fig. 9. — The H-R diagram of the p Oph star forming region outside (upper panel) and within 
(middle panel) the cloud core, as defined by the dashed boundary in Figure ||. Low-mass sources 
observed by WGM and Luhman et al. (1997) (LLR) are also shown. The lower panel compares 
sources in our sample with and without signs of disk activity. The theoretical evolutionary tracks 
of DM97 are given, where the horizontal solid lines are isochrones representing ages of 0.3, 1, 3, 10, 
30, and 100 Myr and the main sequence, from top to bottom. The dashed line in the H-R diagram 
for the cloud core represents Kdereddened = 11, which is the approximate completeness limit of the 
spectroscopic sample for Ak < 3. 
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Fig. 10. — H — K vs. K for the area a) within the p Oph cloud core, as defined by the dashed 
boundary in Figure |l], b) outside the cloud core and within Figure |l|, c) an off field observed by 
BKLT equal in area to the cloud core, and d) fields in the cloud core outlined by dotted lines in 
Figure |^ which were observed at J, H, and K. In a) and b), we show sources classified through our 
K-band spectroscopy. Additional sources have been observed by WGM (K-band spectra) , Luhman 
et al. (1997) (LLR) (optical spectrum), CRCTL {ISO photometry), and BKLT and CRBR [JHK 
photometry). The KLF (entire cloud core) and the HLF (subset of core) are consistent with a 
gaussian reddening distribution {Ak ~ 4.4, a ~ 1.2) towards the background stars behind the 
cloud core, which is simulated in c) by reddening the off-field population. Consequently, the region 
enclosed by the dotted lines is where most background stars are expected to appear in a) and d). 
The dashed lines in d) are the theoretical K magnitudes of substellar objects at 1 Myr (Burrows 
et al. 1997). 
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Fig. 11. — The solid histogram is the IMF derived from the cloud core spectroscopic sample 
with the evolutionary tracks of DM97. All sources with featureless spectra are excluded from the 
IMF. The dashed histogram includes sources that 1) have uncertain late-type spectral types, 2) 
are identified as likely cluster members in ISO observations by CRCTL, or 3) lack spectra and are 
an excess above the background population in the KLF in Figure 0. The complet eness limit is 
indicated by the vertical dashed line and a lower limit to the substellar IMF is provided below this 
boundary. For reference, the field star mass functions of Miller &: Scalo (1979) and Scalo (1986) 
are given. The two lowest and two highest mass bins were given widths of log M = 0.4 due to 
uncertainties in mass estimates. 



-41 - 



n I n 

cloud core: 

BKLT/CRBR photometry 

new photometry 

background model 

spectral types 
uncertain types 
featureless 
WGM 

WGM (background) 
LLR 

CRCTL ISO 



40 - 



30 - 



X5 



3 



20 - 



o 
® 



;1 
•■'l 



10 - 



o 
o 



o 
o 
o 









® 






® 






o 




® 


o 




A 


o 




o 


o 




o 


o 








8 



10 



11 



12 13 
K 



14 



15 



16 



17 



18 



19 



Fig. 12. — The i^-band luminosity function for the cloud core of p Oph. The dotted line is a 
model for the distribution of background stars. 
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Fig. 13. — The pre-stellar clump mass function measured in p Oph from millimetter continuum 
observations (Motte et al. 1998) and the cloud core IMF after a completessness correction assuming 
an age of 0.5 Myr (dashed line in lower panel of Figure 11). The completeness limits of both studies 
are near 0.1 Mq. 
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Fig. 14. — The K-hand spectrum of IRAS64A is compared to data for several FU Ori objects and 
a typical field K giant (reddened by Ay = 20.) Spectra are normalized at 2.2 /xm with constant 
offsets. 
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Table 2. SED Classifications and Spectroscopic Data 



TTV 

ID 


a(2-10 fim)^ 
observed dereddened 


Glass 


tki Br7, a(2-20 ixmy 


GY6 


1.2 






> 2 2 OibO 2 A a20 = 13 


GYlll 


1.2 






1-3 1 8=b0 4 A a20 = 09 


GY214 


1.2 






> 2 3 0±0 1 A a20 = 09 


GY254 


1.0 






> 1, 6.2±0.4 A, a20 = 0.6 


GY265 


0.7 






> 1 1 5±0 7 A a20 = 09 


GY269 


1.4 






1 o20 = 16 


GY274 


1.4 






> 1 2 OzbO 6 A a20 = 11 


GY378 


2.7 






> 2, 4.3±0.4 A, a20 = 2 


GY205 


0.3 




flat 


> 5 5 9±0 4 A a20 < 4 


GY224 


0.2 




flat 


> 5 8 5±1 A a20 = 01 


GY227 


0.2 




flat 


> 0.5 


GSS26 


0.3 


-0.4 


flat/II 


4, 0.75, a20 = -0.1 










2 9±0 2 1 7±0 3 A 


GY21 


0.0 


-0.5 


flat /II 

±±CIjU / J- J- 


1^ 2 0±0 7 A 


GY51 


0.3 


~ -0.4 


flat /II 


2 5 5 OibO 5 1 5±0 5 A 


GY279 


-0.2, 


0.3 -1.0, -0.5 


flat /II 


75 l'^ a20 = 0^ 


GY315 


0.2 


0.9 


flat /II 


1-3 2 OibO 6 A a20 = 01 


GY81 






flat /II? 


1-2 3 OibO 5 A 


EL24 


-0.7 


-0.9 




1-3 2 7±0 2 A a20 = -0 6 


GSS28 


-0.8 


-0.9 


jj 


5 1 9±0 15 A 


GSS29 


-0.8 


-1.2 


jj 


a20 < — 8 


GY23 


-0.5 


-0.9 


jj 


75 a20 = -0 6 


GYllO 


-1.1 


-1.5 


jj 


25 6±0 2A 


GY116 


-0.7 


-1.3 


jj 


75 2 libO 15 A 


GY129 


-1.1 


-1.4 


jj 


> 1 13 7±0 5 A 


GY167 


-0.6 


-0.7 


jj 


1-3 4 7±0 15 A a20 = -0 3 


GY168 


-0.6 


-0.8 




2, 1.3±0.2,a20 = -0.5 


GY211 


-1.0 


-1.4 




0.75, a20 < -0.3 


GY244 


0.3 


_l?g 


II? 


0, a20 < 0.5 


GY247 


-0.9 


-1.5 




0.25, a20 < -0.9 


GY267 


-1.1 


-1.5 




0.25 


GY273 


-0.3 


-1.0 




0, a20 < -0.1 


GY292 


-0.8 


-1.1 




1, 1.0±0.4, 3.0±1 A 


GY308 


-0.7 


-1.1 




0.75, 0.8±0.2 A 


GY314 


-0.8 


-1.1 




0.5, 1.6ib0.3, <0.5 A 


GY319 


-1.3 


-1.3 




0.75, a20 = -1.1 
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Table 2 — Continued 



ID 


o(2-10 iimY 


Class'' 


VK, Br7, a(2-20 /xm)^ 




observed 


dereddened 






nv4nn 

Vjr 1 ^UU 


-1.1 


-1.1 






ST? 4 


-1.0 


-1.1 




1 9 1 ^ (790 — — n 7 

a 7_|_n Q A c_|_n c c c i n c A 


TRS9 


-1.4 


-1.6 




n 


VV kJJ-JUU 


-0.4 


-0.5 




n 7^ 




< -1.9 


< -2.0 









-2.1 


-2.3 




\Jy (XZiU \ ±.\J 




1.5*^ 


1.9^ 


III 


\ J- 


GY12 


< -2.0 


< -2.6 


III 


n 


vjr J. J. 1 


-1.9 


-1.9 


III 


n q_)_n o A 


nY2q 






III?^ 


n 


nY'^7 

V_T 1 O 1 






III?f 


? 


VJ J. Ocf 






III?^ 


? 


OY84 






III?^ 




HYl 


< -2.1 


< -2.6 




n 


HYI 

Vjr 1 -LOO 


-0.9 


-1.7 




n 1 7J_n 4 A 


CYI ^fi 


< -2.2 


< -2.9 


Ill 


n 


GY}72 


< 0.6 


< 0.0 


III? 


n «20 <: 1 1 


nY9'in 

Vjf 1 ^ou 


< -2.10 


< -2.12 




u 


nY9'iS 


< -1.2 


< -2.1 


Ill 


n 


nY9fi9 


-1.5 


-2.3 




1 o_|_n q A 








III?'' 


n 


OY2QB 

V_T J. ZjU^ 


< -2.2 


< -2.4 




n 




< -1.35 


< -1.73 


Ill 


n 


nY'^i n 






III?'' 










III?^ 




(or I 4iU 


< -1.7 


< -2.0 




U 


GY93 


-1.7 


-1.8 







LFAM8 






III?^ 





ROXs39 


-2.6 


-2.6 







SR20 


-1.6 


-1.8 




1 


VSSGll 


< -1.12 


< -1.62 







VSSG22 


< -1.6 


< -2.2 







SR22 


-1.6 


-1.7 




0.25, 1.9±0.4 A 


WSB38 


-1.6 


-1.9 




0.75, 0.6zb0.15 A 


162559-242124 






III?^ 






-46- 



Table 2 — Continued 



ID 



a(2-10 //m) 



a 



Class^ 



rK, Br7, a(2-20 ^m)= 



observed 



dereddened 



162615-241924 
162618-242416 



0.2' 



-2.0 



III?f 
III? 





CRCTL 



'^Data taken from Elias (1978), Lada & Wilking (1984), Young, Lada, & 
Wilking (1986), Wilking, Lada, & Young (1989), and Greene et al. (1994). 

^SED classification (Greene et al. 1994) based on a(2-10 /xm) and the other 



^Continuum veiling and Wx(Br^) measured in the iC-band spectra and the 
spectral index from 2-20 fim when available. 

'^Greene & Lada (1997). 

^Measurement of 20 /xm photometry was at the same epoch as the 10 /um 
data used in alO = 0.3. 

^Tentative SED classification based on JHK colors and veiling in the K- 
band spectra. 

^Dereddening is based on the assumption that H — K \s not affected 
significantly by excess emission, as suggested by rx ~ measured in the 
spectrum. 



data given (see § [4.4.11) . 



^^0(2-6.0 //m). 



a(2-4.5 /xm). 



TABLE 1 

Data for Spectroscopic Sample in p Ophiuchi 



BKLT 


Other ID 


a(1950) 


5(1950) 


previous 


IR. 


.-1 ^-d 
adopt 


^cff 




-^bol 


J — H 


H — K 


K 




Br 






















(K) 














(A 




SR22 


22 


22.7 


22 


55 





K7-M0(GM),M0(CK),M3.5(M) 


K6-M2 


MO 


3850 


0.6 


0.76 


0.93 


0.45 


8.70 


0.25 


1.9±0.4 






22 


35.59 


08 


57 


6 


XV t \j iVl ) 






















162556-242050 


SR4,BL13, YLW25,IRS12 
ROXs6,WSB25,DoAr20 


22 


54.82 


14 


04 


5 


K6(BA),K7(CK) 


K5-A'I2 


K8 


3955 


0.5 


1.9 


1.06 


0. 72 


7.25 


1.5,2.1.5 


6.7zb0.3 
5.5±0.5 


162557-243032 


GSS20,IRS13,ROXs7 


22 


56.11 


23 


46 


6 


K7(GM),K7(BA) 


K5-K7 


K6 


4205 


0.9 


1.1 


1.12 


0.63 


8.37 








162558-242203 




22 


57.23 


15 


17 


7 




K0-K2III: 


b? 








2.57 


1.46 


11.41 


7 





162559-242124 




22 


58.32 


14 


38 


3 




M4 


M4 


3180 


0.0 


0.21 


0.65 


0.45 


9.52 








162601-242521 




23 


00.01 


18 


36 


2 




K-M: 










1.90 


1.23 


11.10 


7 





162603-242337 


GSS23,DoAr21,EL14 
ROXs8,VSS93 


23 


01.65 


16 


52 


4 


KO(BA),KO(M) 


K0-K2 


Kl 


5080 


1.6 


18 


1.26 


0.59 


6.16 








162605-242451(c) 




23 


03.64 


18 


05 


9 


carly(WGM) 




b 








>1.1 


1.91 


13.99 






162609-243411 


GSS25,VSS92,EL16 
YLW28.SR3 


23 


07.84 


27 


26 


6 


AO(SR) 


<F3 


AO 


9520 


2.3 


160 


0.86 


0.49 


6.42 








162610-242056 


GSS26 


23 


08.97 


14 


11 


6 




K5-A'I2 


K8 


3955 


6.0 


1.3 


3.34 


2.25 


9.38 


4,0.75 


2.9±0.2 


162615-241924(c) 




23 


14.49 


12 


40 







M0-M3 


Ml. 5 


3595 


4.6 


O.SO 


2.51 


1.46 


9.98 








162616-242225(c) 


GSS29,EL18 


23 


15.46 


15 


40 


7 


K7-M0(GM:) 


K5-K7 


K6 


4205 


3.0 


3.1 


1.88 


1.04 


8.19 








162617-242023(c} 


GSS28,DoAr24,EL19 
R.OXslOA,WSB27 


23 


15.71 


13 


39 


3 


K7-M0(GM).K5(BA) 


K6-M2 


MO 


3850 


0.8 


1.2 


1.05 


0.58 


8.09 


0.5 


1.9±0.1 


162618-242818 




23 


16.61 


21 


33 


6 


early (WGM) 




b 








2.57 


1.43 


12.93 






162618-242611(c) 




23 


17.54 


19 


27 


4 


M7.5{WGM) 


M5-M6^ 


M5.5 


2925 


3.1 


0.050 


1.75 


1.20 


12.18 






162618-242416(c) 




23 


17.56 


17 


32 


2 


M5(WGM) 


K-M 


M5 


3010 


5.0 


0.097 


2.46 


1.95 


11.94 






162621-242306(c) 


GY6,GSS30.EL21 


23 


19.99 


16 


22 


1 




? 










3.06 


2.51 


8.32 


>2 


2.0±0.2 


162621-242602(c) 


GY5 


23 


20.27 


19 


17 


7 


M7(WGM) 


M5-M6^ 


M5.5 


2925 


1.4 


0.11 


1.13 


0.66 


10.91 


< 1 




162622-242354(c) 


GYIO 


23 


20.78 


17 


10 


2 


M8.5(WGM) 


M6-M7^ 


M6.5 


2780 


4.0 


0.067 


2.16 


1.34 


12.25 






162622-242409(c) 


GYll 


23 


20.89 


17 


25 


1 


M6.5(WGM) 




M6.5 


2780 


1.4 


0.0028 


1.15 


1.22 


14.15 






162622-242254(c) 


GY12 


23 


20.99 


16 


10 


8 


M:2(GM) 


K6-M1 


K8 


3955 


4.6 


0.68 


2.75 


1.87 


9.60 








162623-242 101 (c)S 


GY20A,DoAr24E,GSS31 
EL22,ROXslOB,WSB30 


23 


22.01 


14 


17 


6 


K0(BA} 


K0-K2 


Kl 


5080 


1.8 


10 


1.52 


0.93 


6.44 





0.7±0.2 


162623-242 101 (c)S 


GY20B,DoAr24EB 


23 


22.01 


14 


17 


6 




? 










1.52 


0.93 


6.44 


>2 


1.5±0.4 


162623-244311 


GY17,DoAr25 


23 


22.04 


36 


27 


7 


K4(GM),K7(M) 


K5-K7 


K6 


4205 


0.0 


1.0 


0.93 


0.79 


7.57 





0.9±0.2 


162623-242441 (c) 


GY21,LFAM3 


23 


22.14 


17 


57 


4 


K-M(GM),MO(GL) 


K-M: 


MO 


3850 


3.8 


0.36 


2.45 


1.69 


9.94 


1(GL) 


2.0±0.7 


162624-242449(c) 


GY23,S2,GSS32 
EL23 


23 


22.64 


18 


05 


6 


K(GM) 


K5-M2 


K8 


3955 


3.1 


3.9 


2.15 


1.48 


7.20 


0.75 





162624-241616 


EL24,YLW32 


23 


22.85 


09 


32 







G5-K7 






2.3 




1.95 


1.46 


6.77 


1-3 


2.7±0.2 


162625-242325(c) 


GY31 


23 


23.86 


16 


41 


7 


M5.5(WGM) 




M5.5 


2925 


16 


1.0 




3.98 


13.09 






162625-242446(c) 


GY29 


23 


23.87 


18 


02 


7 




M4-M6 


M5 


3010 


6.3 


0.47 


2.95 


2.11 


10.86 








162627-242643 (c) 


GY37 VSSG29 












]V[6( WGM) 






















162629-241908 


LFAM8 




28 31 


12 


24 






M0-M3 


Ml. 5 


3595 


6.6 


0.46 


3.37 


2.05 


11.19 












23 


28 45 


31 


57 


























±O^OijU- ^ o t^c ) 


^ 1 i- T V 00 'j'^ i 


23 


29 12 


16 


15 






K3 Ml^ 


















r n-Un 
j.0± .5 


16263 1-24253 1(c) 




23 


29.95 


18 


48 





]VI6(^VGM} 


M4-IVI6 


Ms' 


3010 


3 4 


n nqq 


1 86 


121 


11 68 








GY64 


23 


31.11 


19 


53 


4 


iVL ^ W iVL J 




M7 


2720 


3 2 


014 


187 


1 43 


13 33 






162633-243625 


GY65 


23 


31.34 


29 


42 


4 




M2-3III 


b 








2.58 


1.49 


8.91 




Q 


162634-242330(c) 


GY70,S1,GSS35,VSS26 
EL25,ROXsl4,YLW36 


23 


32.78 


16 


47 


1 


B2-B5(BA,GSS,CK) 


<F8 


B3 


18700 






1.56 


0.96 


6.32 








162637-242302 


GY81,VSSG4 


23 


36.42 


16 


19 


4 




K5-M2 


K8 


3955 


4.6 


0.20 


2.79 


1.91 


10.92 


1-2 


3.0±0.5 


162638-242324 


GY84 


23 


37.43 


16 


41 


6 


M6(WGM) 


M0-M3 


Ml. 5 


3595 


4.6 


0.24 


2.40 


1.26 


11.63 








162641-244015 


GY93,WSB37 


23 


39.72 


33 


33 





K-M(GM) 


M4 


M4 


3180 


0.6 


0.31 


0.98 


0.55 


9.20 








162642-242224 


GY107 


23 


41.32 


15 


42 





M3?(WGM) 




b 








2.02 


0.97 


11.40 






162642-242031 


GY110,VSSG2,EL26,GSS37 


23 


41.48 


13 


49 





M0(GM),K6(M) 


K6-A'I2 


MO 


3850 


2.9 


3.5 


1.78 


0.93 


8.00 


0.25 


0.6±0.2 


162643-241635 


VSSG11,IRS19 


23 


42.42 


09 


52 


9 


M3(GM) 


K6-M2 


MO 


3850 


4.0 


1.1 


2.27 


1.27 


9.58 








162644-243447 


GY111,WL12,YLW2 
VSSG30 


23 


42.66 


28 


04 


7 




G5-K7 










3.79 


2.89 


10.18 


1-3 


1.8±0.4 


162645-242309 


GY116,GSS39,EL27 
VSSG28 


23 


43.62 


16 


26 


9 


m:(gm) 


K5-M2 


K8 


3955 


4.7 


2.1 


2.64 


1.62 


8.88 


0.75 


2.1±0.1 


162646-241203 


WSB38,VSS27,ROXsl6 


23 


45.16 


05 


20 


9 


K(GM),G9:(BA) 


K5-M2 


K8 


3955 


1.9 


3.4 


1.47 


0.82 


7.51 


0.75 


0.6±0.1 


162648-242840(c) 


GY128,WL2 


23 


47.05 


21 


57 


9 


K-M(GM) 


K-M: 










>2.95 


3.06 


10.99 


7 





162649-243823 


GY129,WL18 


23 


47.47 


31 


41 


3 




? 






2.6 




1.99 


1.39 


9.82 


>1 


13.7±0. 


162649-242005 


GY135,VSSG3,LFAM18 


23 


47.86 


13 


23 





K7-M0(GM) 


K5-K7 


K6 


4205 


4.4 


4.4 


2.33 


1.20 


8.62 








162651-243242(c) 


GY141 


23 


49.82 


26 


00 


4 


M8.5(LLR) 




M8.5 


2500 


0^- 


0.0026^^ 


0.76 


0.50 


13.87 






162653-242229 


GY150,GSS40,VSSG6 
IRS21 


23 


52.47 


15 


47 


6 




KO-IIII 


b 








2.86 


1.67 


9.85 







162654-242622 


GY153,VSSG5 


23 


53.04 


19 


40 


2 


M;2(GM) 


K6-M2 


MO 


3850 


5.8 


1.2 


3.02 


1.82 


9.88 


0.25 


1.7±0.4 


162654-242231 


GY156,VSSG6,GSS40 


23 


53.58 


15 


49 


4 




K6-M1 


KS 


3955 


6.1 


1.3 


3.07 


1.78 


10.19 








162656-242838(c) 


GY163 


23 


55.47 


21 


56 


8 


M2.5(WGM) 




M2.5 


3430 


9.0 


0.26 


>2.06 


2.28 


12.66 






162658-244529 


GY168,SR24N,IRS24 
EL28 


23 


56.58 


38 


48 


1 


K2(GM),M0.5(CK),M0(M) 


K5-A'I2 


K8 


3955 


2.0 


2.3 


1.70 


1.18 


7.42 


2 


1.3±0.2 


162658-244534 


GY167,SR24S,IRS25 
EL28 


23 


56.64 


38 


52 


9 


K2(GM),K2:(CK),K7(M) 


G5-K7 


K2 


4900 


1.2 


2.4 


1.50 


1.17 


7.08 


1-3 


4.7±0.1 


162659-243556(c) 


GY172,WL14 


23 


57.42 


29 


14 


7 




M4 


M4 


3180 


5.1 


0.14 


2.74 


1.64 


11.74 
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TABLE 1— Continued 



BKLT^ 


Other ID^ 


c.(1950) 


5(1950) 


previous*' 


IR 


adopt*^ 


T ff 
(K) 


A J 


-^bol 

f 7",/-N 1 

K^Q > 


J — H 


H - K 


K 




Br-7 

(A) 


1 62702-243726 (c) 


G Y182. WL16 


24 


00.64 


30 


45.0 




B8- A7^ 










3.56 


2.66 


7,92 


7 


7.2±0.2 


1 62704-242830 (c) 


GY192.WL1.YLW6 


24 


02.61 


21 


49.3 




>M3 










3.15 


2.26 


10.81 


7 







IRS29.VSSG31 
































16 270 6 -242837(c) 


GY202 


24 


04.47 


21 


56.0 


Jvl7( W CjJvI ) 


JVlD-iVl ( 


M6.5 


2780 


3.9 


0.022 


2.09 


1.70 


12.97 






162 706- 2438 14(c) 


G Y 2 5 . W L 1 7 


24 


05.06 


31 


33.3 




I' 










3.33 


3.29 


10.28 


>0.5 




162709- 243408(c) 


G Y2 1 1 . WLIO 


24 


07.48 


27 


27.0 




K5-M2 


K8 


3955 


3.6 


1.4 


2.22 


1.39 


8.85 


0.75 


1 .4±0.2, < 0.8 


162709- 2437 18(c) 


G Y2 14,EL29,^^L15 


24 


07.71 


30 


37.4 














5.20 


4.47 


7.54 


>2 


3.0±0.1 




































162710-241914 


SR21, VSSG23,VSS31 


24 


08.78 


12 


33.7 


G1(M) 


P0-F8 


F4 


6590 


2.7 


51 


1.23 


1 .03 


6,30 


<; 1 


Q 




EL30, YLW8 
































162711-244046 


GY224 


24 


09.38 


34 


05.7 




? 










>3.31 


2.90 


10.79 


>0.5 


8.5±1.0 


162711-243S32(c) 


GY227 WL19 


24 


10.00 


31 


51.2 




7 












>4.27 


11.23 


>0.5 


Q 


162711-242343 


GY228 IRS32 


24 


10.25 


17 


02.6 














2 63 






7 




162 7 12- 243449(c) 




24 


10.49 


28 


08.4 




K-M: 










2.37 


1. 76 


11.53 


7 


Q 


162713-244133 


G Y232 


24 


11.44 


34 


52.5 






]-, 


















162714-245132 


WSB45,R.OX.s20A 


24 


12.34 


44 


51.7 


M 5 ( B A ) 


M1-M4 


PMS?.dMc? 








0.79 


0.34 


10.39 


Q 


Q 


162715-245137 


WSB46,R.OXs2nB 


24 


12.98 


44 


57.2 


M 2 ( B A ) 


M0-M3 


PMS? dMc? 








0.90 


0.46 


9.51 








162715-242640 


GY239 inS34 YLW^lOB 


24 


13.99 


19 


59.9 


K ( G M ) 


















7 




16271 5- 243843 (c)S 


GY240A,WL20 


24 


14.08 


32 


02.9 


K7-M0(GM) 


K6-M2 


MO 


3850 


4.9 


1.4 


2.67 


1.57 


9.21 


0.25 





16271 5- 243843 (c)S 


GY240B, WL20 


24 


14.08 


32 


02.9 


K-M(GM) 


K5-K7 


K6 


4205 


4.9 


0.55 


2.67 


1.57 


9.21 








16271 7-242856 (c) 


GY244,YIjW^12 A,IR,S37 


24 


16.06 


22 


16.1 




M4 


M4 




--^ llJ 




>3.1 


2.95 


10.95 








V bb(j^D 
































162718- 242853(c) 


GY246,WL5,VSSG26 


24 


16.65 


22 


12.6 


F7(GM) 


<G0 










>2.68 


4.04 


10.28 


7 







IRS38 






























1627 18- 242906(c) 


Lr Y 247 , W lj4,iRbdy 


24 


16.95 


22 


25.6 




M0-M3 


Ml. 5 


3595 


5.4 


1.6 


3.00 


1.93 


9.13 


0.25 




V Oblj^O 






























162719-244 139(c) 


/-I veycrt o'Dio T'DC^in 


24 


17.67 


34 


59.5 




K6-M2 


MO 


3850 


0.2 


0.66 


0.73 


0.28 


8.41 






RUA-S^I , Y LW 13A 










JVll l^OlS. ) 






















162721-244 142(c) 


Ij Y 202, I JjW Id ±3 ,i±t,b4^ 


24 


19.63 


35 


02.5 














3.90 


2.90 


8.41 


> 1 




162721-244335 


GY253 


24 


19.95 


36 


55.5 




M2-M4 


M3 


3350 


8.2 


0.83 


4.03 


2.53 


10.78 








16272 1-242953 (c) 


GY254, Y]JW14,W^L6 


24 


20.26 


23 


12.9 


■ ■ ■ 












>2.61 


4.35 


10.04 


>1 


6.2±0.4 


162722-241759 


VSSG22 


24 


21.51 


11 


19.1 


K7-M0(GM) 


K5-K7 


K6 


4205 


4.6 


1.8 


2.48 


1.38 


9.41 








162726-243923 (c) 


GY262 


24 


24.67 


32 


44.0 




K7-M3 


Ml 


3680 


6.4 


1.0 


3.46 


2.30 


9.77 


^ 


2.2±0.3 


162726-244051 (c) 


G Y265 , YLWl 5 A , IRS43 


24 


25.09 


34 


11.6 














>3.83 


3.71 


9.46 


> 1 


1.5±0.7,<0.5 


162727-243116 


GY267,VSSG25,EL31 


24 


25.83 


24 


37.0 


K-M(GM) 


K6-M2 


MO 


3850 


3.1 


0.92 


1.93 


1.11 


9.30 


0.25 





162728-243934 


WL13 

Ij Y 20ti, I JjW lDA,iH,b44 


24 


26.22 


32 


54.5 






t5 








>3.91 


3.44 


9.65 






162728-242721 


GY273,VSSG18,EL32,IRS45 


24 


26.92 


20 


42.1 




K5-M0 


K6.5 


4132 


5.0 


0.86 


2.99 


2.09 


9.39 








162729-243917 


GY274, Y1JW16B ,IRS46 


24 


27.63 


32 


37.2 














>2.35 


3.19 


11.46 


> 1 


2.0±0.6 


162730-242744 


GY279,VSSG17,EL33 


24 


28.65 


21 


04.6 


M (GM) , > M2(GL) 


K5-M2 


K8 


3955 


6.9 


1.8 


3.80 


2.69 


8.95 


0.75,1(GL) 







IRS47 
































162730-242457 


GY284 


24 


29.37 


18 


17.4 




M0-M3 


Ml. 5 


3595 


2.4 


0.35 


1.69 


0.97 


10.04 








162733-2441 15(c) 


GY292 


24 


31.29 


34 


36.2 




K5-M2 


K8 


3955 


2.9 


2.1 


2.04 


1.37 


7.92 




1.0±0.4,3.0±1 


162735-243834 


GY295 


24 


33.50 


31 


55.0 




M4 


M4 


3180 


1.3 


0.36 


1. 13 


0.51 


9.64 








162736-242833 


GY297 


24 


34.95 


21 


54.7 




M0-M3 










0.68 


0.25 


11.04 








162737-243035 


GY304, YL,W46 A,IRS48 












<F3 










1.88 


1.23 


7.42 







162738-243043 


GY306,IRS50 


24 


36.57 


24 


03.9 




M4 


M4 


3180 


3.1 


0.63 


1.87 


1.01 


9.59 








162738-243658 


G Y308 , IRS49 


24 


36.67 


30 


19.1 




K5-M2 


K8 


3955 


2.9 


1.7" 


1.96 




8.31 


0.75 


0.8±0.2,<2 


162738-243838 


GY310 


24 


36.90 


31 


59.1 


M8.5(WGM) 




M7 


2720 


1.7 


0.079 


1.29 


0.83 


11.08 






162739-243914 


GY314 


24 


37.68 


32 


35.6 


M(GM) 


K6-M2 


MO 


3850 


2.0 


1.4 


1.55 


0.92 


8.35 


0.5 


1.6±0.3,<0.5 


162739-244316 


GY315,YLW45,IR,S51 


24 


37.97 


36 


37.0 


■ ■ ■ 


G5-K7 






9.5 




4.70 


3.49 


8.93 


1-3 


2.0±0.6 


162740-242205 


GY319.SR.9.EL34.R,OXs29 


24 


38.83 


15 


26.4 


K4(GM),K7(CK!) 


K5-M2 


K8 


3955 


0.0 


2.0 


0.77 


0.51 


7.20 


0.75 







IR,S52 , Do Ar34. WSB54 










K4/K6(B A) 






















1 62742-243849 




24 


40.96 


32 


10.9 


]VI4 ( WGM) 




M4 


3180 






1 78 


88 


10 59 






162746-242323 


GY351.VSSG13,EL35 


24 


45.25 


16 


45.2 




MOIII 


b 








2.20 


1.12 


7.29 




Q 


162749-242540 


GY372,VSSG14,EL36 


24 


48.37 


19 


02.5 


A7(GM),A7:(M) 


<F8 










1.20 


0.82 


7.32 








162751-243145 


GY378,YLW52,IRS54 


24 


50.26 


25 


07.6 




? 










3.13 


2.63 


10.87 


>2 


4.3±0.4 


162755-242619 


GY400,SR10 


24 


53.97 


19 


41.0 


M1.5(CK} 


K6-M2 


MO 


3850 


0.0 


0.55 


0.50 


0.60 


8.74 


0.5 


2.5±0.3 


162757-244004 


GY410 


24 


56.08 


33 


26.1 


M6(GM) 


K6-M0 


K7 


4060 


2.7 


0.64 


1.72 


0.92 


9.78 








162757-243751 


GY411 


24 


56.19 


31 


13.3 




G8.5III 


b 








2.20 


1.21 


9.52 







162803-242632 


GY453,VSSG16,EL37 


25 


02.13 


19 


55.5 




K2III 


b 








1.10 


0.77 


6.89 







162809-242322 


IRS 58 

EL38,VSSG15 


25 


07.71 


16 


45.0 




K2-3III 


b 








1.52 


0.82 


7,12 







162816-243657 


WSB60,YLW58 


25 


14.88 


30 


20.9 


K-M(GM) 


K5-M2: 


K8 


3955 


0.7 


0.28 


1.16 


0.78 


9.43 


0.75 





162823-242241 


SR20,ROXs33,DoAr38 


25 


21.82 


16 


05.7 


GO:(BA) 


K5-M2 


K8 


3955 


2.1 


0.93 


1.69 


1.13 


8.55 


1 







WSB61 
































162923-241359 


IRAS64A 


26 


22.00 


07 


27.1 














1.79 


0.91 


6.45 


7 







ROXs39 


27 


33.8 


27 


51.0 


K5(BA) 


K5-K7 


K6 


4205 


0.0 


1.4 


0.66 


0.27 


8.00 









TABLE 1— Continued 



BKLT* 


Other id'' 


a(1950) 
min sec 


5(1950) 
/ // 


previous'' 


IR 


adopt'* 


Tetf 
(K) 




^bol 


J - H 


H - K 


K 




Br^ 
(A) 





^Sources falling within the cloud core as defined in Figure 1 are denoted with (c). 

^ SR=Struve &; Rudkj0bing (1949), DoAr=Dolidze & Arakeylan (1959), S and GSS=Grasdalen, Strom, & Strom (1973), VSSG=Vrba et al. (1975), 
VSS=Vrba, Strom, &c Strom (1976), EL=Elias (1978), ROX=Montmerle et al. (1983), WL=Wilking &c Lada (1983), YLW=Young, Lada, &: Wilking 
(1986), WSB=Wilking, Schwartz, &c Blackwell (1987), IRAS=Table 1 of Wilking, Lada, &c Young (1989), IRS=Table 2 of Wilking, Lada, & Young (1989), 
LFAM=Leous et al. (1991), GY=Greene & Young (1992). 

*^IR spectral types of Greene Meyer (1995), Greene & Lada (1997), and Wilking, Greene, & Meyer (1998). Optical spectral types of Struve & 
Rudkj0bing (1949), Grasdalen, Strom, & Strom (1973), Cohen &: Kuhi (1979), Bouvier &c Appenzeller (1992), Luhman, Liebert, &c Rieke (1997), and 
Martin ct al. (1998). 

'^Background and foreground stars are indicated by "b" and "f ' . 

^Derived from the IR spectra in this work. 

^Data from Williams et al. (1995) used in classification. 

®GY20 and GY240 are each binary systems which appear as one source in the coordinates and photometry of BKLT. 

^Luhman, Liebert, &: Rieke (1997). 

'JT-band classification by Biscaya et al. (1999). 

■•Extinction was derived by assuming that — if is not affected signicantly by excess emission, as suggested by the rj^ measured in the spectrum. 
^The optical spectral type is M8-M9III while the IR spectrum exhibits CO absorption equivalent to K4III and is very similar to spectra of FU Ori 
objects. 

NOXE. — JHK photometry and coordinates are taken from BKLT. The data for ROXs39 and SR22 are from Greene et al. (1994) and references therein. 
Photometry for GY31 from Comeron et al. (1993). Coordinates are prefixed by 16^ and —24°. 



